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Abstract. We study soft hadron production in relativistic heavy ion collisions in a wide range of reaction 
energy, 4.8 GeV < t/snn < 200 GeV, and make predictions about yields of particles using the statistical 
hadronization model. In fits to experimental data, we obtain both the statistical parameters as well as 
physical properties of the hadron source. We identify the properties of the fireball at the critical energy 
threshold, 6.26 GeV < -J s^ N < 7.61 GeV, delineating for higher energies hadronization of an entropy rich 
phase. In terms of the chemical composition, one sees a phase which at low energy is chemically under- 
saturated, and which turns into a chemically over-saturated state persisting up to the maximum accessible 
energy. Assuming that there is no change in physical mechanisms in the energy range 15 > -^/snn > 
200 GeV, we use continuity of particle yields and statistical parameters to predict the hadron production 
at y^snn = 62.4 GeV, and obtain total yields of hadrons at ^/snn = 130 GeV. We consider, in depth, the 
pattern we uncover within the hadronization condition, and discuss possible mechanisms associated with 
the identified rapid change in system properties at ^/sJf N . We propose that the chemically over-saturated 
2+1 flavor hadron matter system undergoes a 1st order phase transition. 
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1 Introduction 

It is believed that the deconfined phase of matter is 
formed at sufficiently high energy and reaction volume 
reached in the most central collisions of heavy ions at the 
top RHIC energy [1]. The question is where this criti- 
cal energy threshold V^nn is. We pursue this point in 
this systematic study, in order to explore possible phase 
changes occurring as function of collision energy [2]. We 
furthermore compare our results qualitatively to the be- 
havior seen as function of the reaction volume [3] . 

The tool, used in our study of soft hadron produc- 
tion, is the generalized statistical hadronization model 
(SHM) which allows for particle yields to be in full chem- 
ical non-equilibrium [4]. SHM is capable to describe, in 
detail, hadron abundances and has been considerably re- 
fined in past decade, after its formulation by Fermi and 
Hagedorn [5]. 

We present and/or extend here results of analysis of 
the energy dependence of total hadron production yields 
for: 

a) fixed target symmetric Au-Au reactions at the top 
available AGS projectile energy 11.6 AGeV (energy per 
colliding nucleon pair ^/snn = 4.84 GeV), 

b) fixed target symmetric Pb-Pb reactions at SPS at 20, 
30, 40, 80 and 158 AGeV projectile energy. This we refer 



to as SPS energy range, 6.26 < ^/snn < 17.27 GeV, 
c) the Au-Au reactions in the collider mode at RHIC 
in 65+65, 100+100, and also at 31.2+31.2 AGeV reac- 
tions for both total, and central rapidity yields. This is 
the RHIC energy range, 62.4 < ^siiT < 200 GeV. 

Experimental data analysis at RHIC was carried out 
for the central rapidity region yields at ^/snn = 130 and 
200 GeV, and for the full hadron yields at 200 GeV. The 
results we present for total hadron yields at ^/snn = 130 
GeV and 62.4 GeV, arise from our model considerations 
alone. This is also in part the case for for the central ra- 
pidity yields at ^/snn = 62.4 GeV . 

As a first step, we aim to describe at each reaction en- 
ergy the hadron yield data. We obtain in this process the 
statistical hadronization model (SHM) parameters, which 
allow to evaluate the yields of all (also of unobserved) 
particles. One can see SHM analysis as a method of how 
the known experimental hadron yield data can be extrap- 
olated to obtain the unobserved hadron yields. For this 
reason, we also attempt to extrapolate to reaction ener- 
gies and phase space coverage which is not allowing, for 
lack of data, a SHM fit. For example, for the 31.2+31.2 
AGeV case, we interpolate strange particle yields, which 
are known below and above this energy, and/or fix certain 
SHM parameters which show continuity as function of re- 
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action energy, respecting in the process the constraints of 
the SHM. 

In our analysis, we are seeking consistency in the re- 
sults across the reaction energy. This is of importance 
when the number of measurements is not much greater 
than is the number of observables. When we are able to fix 
the values of statistical parameters by consistency check 
across energy range, then the number of used parame- 
ters is reduced and the statistical significance shown in 
this work is for the number of parameters actually fitted. 
However, more often, in the study of statistical signifi- 
cance, we do not account explicitly for consistency across 
energy range. For this reason, in most cases, the statistical 
significance we present is a lower limit. 

Once a statistically significant description of the data 
sample at an energy is achieved, we have available the 
yields of all soft hadronic particles and their resonances. 
We sum partial contributions of each particle species to 
quantities such as entropy, strangeness, baryon number, to 
obtain the properties of the fireball at the time of particle 
production (hadronization). In this way, we evaluate fire- 
ball breakup pressure, entropy, baryon number, strangen- 
ess yield and the thermal energy content. In this approach, 
the kinetic energy content associated with the collective 
flow of matter is not considered — this requires a study of 
particle m± and rapidity spectra, beyond the integrated 
hadron yields. 

High strangeness [6], and entropy content [7,8] of a 
dense hadronic matter fireball are the anticipated charac- 
teristic property of the color deconfined state of matter. 
Once formed, this enhancement of strangeness and en- 
tropy is also the property of the final hadronic state: first 
principles require that entropy must increase in the fireball 
expansion, as well as in the ensuing hadronization process; 
model studies show that once strangeness is produced, it 
remains present during expansion of dense matter, it can 
slightly increase during hadronization [9,10]. 

Particle yields, and pion yield in particular, provide 
natural measure of entropy yield, while the kaon yields, 
and in particular the K+ yield, are an approximate mea- 
sure of the total strangeness yield across all reaction ener- 
gies [7]. The yield ratio K + /n + has been studied as func- 
tion of reaction energy in the SPS energy domain and a 
strong 'horn' like feature has been discovered [2]. 

This suggests a change in the reaction mechanism of 
particle production, occurring in central collisions of Pb- 
Pb, in the energy interval 6.26 GeV < ^/snn < 7.61 GeV, 
the two limits correspond to 20 and respectively, 30 AGeV 
Pb beams on fixed target. This energy range is just at the 
predicted threshold of quark-gluon formation arising con- 
sidering balance of energy deposition and relativistic re- 
action dynamics [11]. Possibility of a rich phase structure 
of the deconfined phase at high baryochemical potential 
and finite temperature further enhances our interest in the 
study of this reaction energy domain [12]. 

To describe experimental results indicating the pres- 
ence of a critical ('cr') energy threshold, one can, in first 
instance, use two different reaction models which apply be- 
low and, respectively, above a postulated energy threshold 



for a phase transformation [13]. However, this presupposes 
the most important outcome, namely that there is an en- 
ergy dependent change in dense matter fireball structure 
at its breakup. Moreover, such an approach does not pro- 
duce as result of analysis an insight into the structural 
change that occurs, and which could be compared with 
predictions. Instead, the structural change is part of the 
hypothesis under which the analysis is carried out. For this 
reason, the methodology we choose here is more general. 

We use in this work the software package SHARE 
(statistical hadronization with resonances) [14], the pub- 
lic SHM suit of programs, where the methods of SHM 
analysis are described in greater detail. Of particular im- 
portance here is that the full mass spectrum of hadron 
resonances is included [15]. SHARE implements two fea- 
tures important for the full understanding of the K+/-7T+ 
horn: 

1) the isospin asymmetry driven by proton-neutron asym- 
metry, which is particularly relevant at low reaction ener- 
gies, 

2) the chemical non-equilibrium (phase space under-satu- 
ration and over-saturation) for strange and light quarks. 
These two features appear to be essential to obtain a de- 
scription of the K + /7r + energy dependent yield. 

We first describe, in next section 2, features of the 
data sample we use, discuss the input data and results of 
the fits for the AGS/SPS and RHIC energy range sepa- 
rately. We discuss the resulting statistical parameters and 
the confidence level of our fits. We survey, in both tabu- 
lar and graphic form, the energy dependence of particle 
yields of interest, including an explanation of the K+/-7T+ 
horn. We then discuss physical properties of the fireball 
at point of chemical freeze-out in section 3, show the en- 
ergy dependence of the model parameters and of physical 
properties, and address the strangeness and entropy pro- 
duction. We discuss the results of our analysis and present 
their interpretation in the final section 4. 



2 Fit procedure and hadron multiplicities 

2.1 General remarks 

The measured experimental results are available for either 
total particle yields, N^, or for central rapidity yields, 
dN/dy. At RHIC energy scale, we will study both data 
sets, though is rather incomplete at this time. At 
AGS and SPS, we will solely consider iV 47r , in order to 
minimize the impact of the shape of the longitudinal un- 
stopped matter flow on the outcome of the analysis. 

At SPS, a semi-distinct central rapidity domain is only 
present in the top SPS case, its re-analysis will make good 
sense once the RHIC dN / dy data extend to the minimum 
accessible energy domain which is close to top SPS energy 
range. However, this will require introduction of models 
of collective matter flow, a step which we do not wish to 
take in this work. At high RHIC energies, we presume that 
the fragmentation regions are sufficiently separated from 
the central rapidity domain as to allow the study of the 
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rapidity particle distributions dN/dy, at central rapidity, 
in a model independent fashion. 

We include, in our consideration of the total particle 
yields N 47T , the trigger condition which defines the par- 
ticipant 'wounded' nucleon number N w . This has to be 
equal to the total net baryon number b = B — B contained 
within the final state particle multiplicities. Furthermore, 
both for N in and central rapidity yields dN/dy, we con- 
sider two constraints: 

a) the fraction of protons among all nucleons (0.39 for 
heavy nuclei) establishes a fixed final ratio of all electrical 
charge Q to the total final state baryon number b - the 
ratio Q/b is preserved in any fraction of a volume of cen- 
trally colliding nuclei; it is a measured quantity given the 
fact that both target and projectile are known; 

b) strangeness (valance s-quarks) content of hadrons prior 
to weak decays has to be (up to systematic experimen- 
tal error) balanced by antistrangeness (valance s-quarks) 
bound in hadrons for the study, and nearly balanced 
when considering the central rapidity dN/dy distributions. 

As our prior studies showed [16], any deviation from 
strangeness conservation as function of rapidity is, in gen- 
eral, smaller than the typical 10% systematic error of the 
experimental data points . It is the level of systematic er- 
ror in the particle yields which determines the precision at 
which we have to assure strangeness conservation. Forc- 
ing exact balance can create an aberration of the fit, since 
the sharp constraint is inconsistent with several indepen- 
dent measurements which contribute to the cancellation. 
For example, at several SPS energies the systematic errors 
between K + and A which control the yields of s and, re- 
spectively, s quarks, do not cancel to better than 8% level. 
This can be checked without a fit in a qualitative study of 
the key particle yields. 

Another reason to be cautious with the strangeness 
conservation is that the spectra of hadrons we are using 
could contain wrong entries (e.g., pentaquark states which 
we in view of recent experimental results do not anymore 
include in the input data set, or wrong spin-isospin assign- 
ments for little known states). Moreover, we maybe miss- 
ing some relevant undiscovered resonances. These effects 
are largest when the baryon asymmetry is largest, since 
the strangeness balance condition probes at large bary- 
ochemical potential the mass spectrum of strange baryons 
and mesons separately, with mesons dominating in an- 
tistrangeness and baryons important in the strangeness 
count. 

For this reason, our strangeness conservation proce- 
dure is as follows: when a first fit shows strangeness a 
slight strangeness asymmetry, we find the best param- 
eters for the fit with a loose, systematic error related 
strangeness conservation constraint allowing, e.g., a 10% 
deviation from balance as a fit input that is we request 
(s — s)/(s + s) = ± 0.1. Since we present confidence level 
and profiles of the fit, and we wish to have from energy 
to energy comparable results, we redo the fit with a fixed 
preferred value of the strangeness fugacity A s as is done 
in case of using exact strangeness conservation. 



In this way, we obtain a data fit with the same mecha- 
nism of approach as for the cases where exact strangeness 
conservation is used to fix one parameter, so that confi- 
dence levels are comparable. We find —0.07 > (s — s)/ (s + 
s) > —0.1 in the SPS energy domain. The asymmetry fa- 
vors an over-count of s quarks in emitted hadrons. It is 
moderate in its relative magnitude, staying within the sys- 
tematic errors of the measurements used in this study. We 
will state the strangeness balance explicitly when present- 
ing the computed particle yields. Note that addition of 
pentaquarks to the hadron spectrum decreases this asym- 
metry by 0-3%, but has otherwise minimal influence on 
the fit results presented. 

As the above discussion of strangeness conservation 
shows, conserved quark quantum numbers introduce yield 
constraints, which particle multiplicities cannot deviate 
from. How a subset of SHM parameters determine a set 
of particle ratios has been shown for the first time in 
1982 [17]. An nice example is the chemical relation be- 
tween the K~/K + and p/p demonstrated experimentally 
in 2003, see figure 4 in [18] , a development based on the re- 
discovery of the SHM constraints in 2000 [19]. Since SHM 
with its chemical consistencies has been very successful 
in helping understand hadron production, we embark on 
further data verifications at each energy, checking the con- 
sistency of experimental data with SHM. 

A suspect particle yield can be further cross checked 
studying the behavior of this particle yield as function of 
energy. Such consideration is very important since we are 
searching for a change in the physical properties of the 
fireball as function of energy, and we do not want the out- 
come to be even in part the result of a statistical fluctua- 
tion in the reported yield of a subset of particles. We find 
inconsistencies (see next paragraph) in the particle yield 
effects. None of these influence decisively the findings we 
report here, in part because of the more lax attitude we 
take toward the constraint on strangeness conservation we 
described above. Moreover, considering the large number 
of experimental data considered, fluctuations in experi- 
mental data sample must occur. 

Specifically, we did not use the A(1520) nor fl and fl 
yields obtained at 158 A GeV in our fit. The preliminary 
yl(1520) value at top SPS is ,4(1520) = 1.45 ± 0.4 [20]. 
This is within 3 s.d. of the SHM yield. However, this ex- 
ceptionally narrow resonance may be subject to additional 
effects [21] and we felt that it is more prudent to not in- 
clude its study here. The experimental yields of fl and fl 
at 158 A GeV are contrary to the 71(1520) larger than the 
SHM model predicts. These particles are produced very 
rarely and for this reason any novel mechanism of produc- 
tion [23] would be first visible in their yield. We believe 
that it is also prudent to not include these in the study, 
even if the deviation from fit would be at 2 s.d. level. 

2.2 AGS and SPS energy range fit 

To assure the reproducibility of our analysis, we will de- 
scribe in detail the input particle yields that are used, 
for AGS/SPS energy domain, and for RHIC domain in 
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the next subsection. The set of particles available at AGS 
arises from several experiments, we have previously re- 
ported in detail the SHM analysis at the top AGS en- 
ergy [24], which input and fit results are restated here. 

The study of AGS results was performed in [24] for sev- 
eral possible cases, such as with and without <f>/K yield, 
strangeness (non)conservation. The results here presented 
are for the case in which the <fi/K yield is fitted and strang- 
eness is conserved. Differences in theoretical fit detail yield 
similar fit result which show the robustness of the ap- 
proach. 

For this work, the analysis of the particle yields 
of the NA49 experimental group available at 20, 30, 40, 
80 and 158 A GeV [25] has been carried out. This work 
extends significantly our prior study of the 40, 80 and 158 
A GeV NA49 done when many fewer measurements were 
available [52]. Moreover, the SHARE package used offers 
additional theoretical features which were not fully im- 
plemented earlier: the consistent description of the yields 
of different charges hadrons (e.g., protons and neutrons, 
7r + and 7r~ by means of A73 allows to fix the net charge 
fraction Q/b. The most relevant difference to the earlier 
study is, however, that we can address the two newly mea- 
sured reaction energies, 20, 30 A GeV. This, along with 
the AGS 11.6 A GeV data, including the recently pub- 
lished 0-yield [53], allows to recognize a major change in 
the behavior of the hadronizing fireball [2] . 

The input data we considered for the AGS and SPS are 
presented in top part of the table 1 . The statistical param- 
eters are seen below these input data. In carrying out the 
data analysis, we use the full grand-canonical statistical 
set of seven parameters: volume V, freeze-out tempera- 
ture T, chemical quark fugacities A 9iS , quark occupancy 
parameters 7 9 and 7 S , and third component of the isospin 
fugacity A13. The fitted values of these 7 parameters are 
seen near bottom of the table 1, which is followed by en- 
tries for the central values of the two chemical potentials: 

^ B = 3TlnA 9 , (1) 
Ms = Mb/3 -Tin A s . (2) 

The uncertainties in the value of statistical parameters 
comprise the propagation of experimental measurement 
error through the fit, as well as ambiguity due to statisti- 
cal parameter correlations arising. In some instances this 
effect is very small, in others rather large. This wide dis- 
parity is possible, as sometimes the data set is sufficiently 
constraining, and in others it is not. The most interesting 
result, we notice in table 1, is the sudden shift in the val- 
ues of the phase space occupancies j q and 7 S observed as 
reaction energy rises from 20 to 30 A GeV. The value of 
chemical freeze-out temperature T changes accordingly to 
counterbalance the effect of a rapid change in "f q and j s 
on some particle multiplicities. We will discuss this change 
in behavior in great detail in what follows. The steady de- 
crease of baryochemical potential /xb with reaction energy 
follows the enhancement in global yield of hadrons. At cen- 
tral rapidity the steady increase of baryon transparency 
with increasing collision energy yields a smaller value of 
fiB- The total particle yields we consider here yield an 



average over the entire rapidity range of /ib • The associ- 
ated value of fis is controlled by strangeness conservation 
condition, as discussed. 

As seen in table 1, we occasionally fix the value of j q . 
The value we choose is the the best value which emerges 
from study of x 2 profile, see figure 1. We fix the best j g 
in order to reduce the correlations between parameters, 
given the small number of degrees of freedom. Exclud- 
ing from the count of parameters \ s which is related to 
(near) strangeness conservation, there are 6 (or 5 when 
7 g is fixed) parameters while there are between 9 and 12 
data inputs at each energy considered in table 1. We thus 
have 3-6 degrees of freedom (d.o.f ) for the fits carried out 
at AGS and SPS. Not all of the NA49 SPS energy range 
results we use are published in final form. 

We show, in figure 1, the reliability of the fits we ob- 
tained at different reaction energies as function of 7 ? , the 
light quark phase space occupancy. The results for AGS 
and SPS are accompanied by those for central rapidity 
RHIC fits we will address below. The top frame, in fig- 
ure 1, shows x 2 /dof . The associated significance level P[%] 
is seen in the bottom frame. We include P[%] as result, 
since the number of degrees of freedom in each fit is small 
and it is hard to judge the significance of a small value of 
X 2 /doi. 

We study dependence of x 2 /dof and P[%] on 7 g since 
we see, in table 1, that the two parameters which un- 
dergo a rapid change as function of reaction energy are 7 g , 
and to a lesser degree, the freeze-out temperature T. The 
rapid change with "f q is prominent in figure 1 where P[%] 
peaks for the lowest two energies (11.6 and 20 AGeV) at 
7 9 < 0.5, while for all other collision energies it grows to 
maximum value near "f q ~ 1.6, where the Bosc singularity 
of the pion momentum distribution 7^ ~ e m i/' r _ 

The reader can see, in figure 1, that setting the value 
7 9 = 1 will yield a set of energy dependent individual fits 
which appear to have a good confidence level. However, 
the energy dependence of the particle yields derived at 
this fixed j q = 1 condition is less convincing. It is the 
rapid shift in the best j q as function of reaction energy 
which allows to describe the 'horn' feature in the K + /tt + 
data (see below figure 2). Without variable 7, this horn 
feature is largely erased, see, e.g., figure 4 in [54], and 
the dashed and dotted lines in figure 2. We will return to 
describe this effect in section 2.4, and discuss this physics 
in more detail in section 4.1. We believe that, in the study 
of energy dependent particle yields, the use of the highest 
confidence level SHM results with 7 g ^ 1 is required for 
the description of energy dependent particle yield data. 

Regarding the weak decay (WD) contributions: in the 
fits to experimental data, we have followed the corrections 
applied to the data by the experimental group(s). For 20 
and 30 GeV in A and A, the data we use includes the WD 
of S, E, Q and fi, these matter mainly in consideration of 
antihyperon yields. At all higher SPS energies all hyperon 
WD decays are corrected for by the NA49 collaboration, 
the decays are always corrected. Similarly, decays of 
kaons into pions are corrected for at all SPS energies. At 
AGS 11 GeV, all yields we consider are without WD con- 
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Table 1. The input total particle multiplicities data at top, and, below, the resulting statistical parameters for AGS and 
SPS energy range. At bottom, we state the chemical potential corresponding to the central values of the fugacities. For each 
projectile energy E [A GeV], we also present in the header the invariant center of momentum energy per nucleon pair, ^/snn 
[GeV], the center of momentum rapidity and the centrality of the reaction considered. The \ s values marked with a * are result 
of a strangeness conservation constraint. 



E[AGoV] 


11.6 


20 


30 


40 


80 


158 


V s nn [GeV] 


4.84 


6.26 


7.61 


8.76 


12.32 


17.27 


Vcm 


1.6 


1.88 


2.08 


2.22 


2.57 


2.91 


N4 n centrality 


most central 


7% 


7% 


7% 


7% 


5% 


R = p/n + , N w 


R = 1.23 ±0.13 


349±6 


349±6 


349±6 


349±6 


362±6 


Q/b 


0.39±0.02 


0.394±0.02 


0.394±0.02 


0.394±0.02 


0.394±0.02 


0.39±0.02 




133.7±9.9 


184.5±13.6 


239±17.7 


293±18 


446±27 


619±48 


R — 1T~ /7T + , 7T~ 


R = 1.23 ±0.07 


217.5±15.6 


275±19.7 


322±19 


474±28 


639±48 


R = K+/K", K+ 


R= 5.23 ±0.5 


40±2.8 


55.3±4.4 


59.1±4.9 


76.9±6 


103±10 


K" 


3.76±0.47 


10.4±0.62 


16.1±1 


19.2±1.5 


32.4±2.2 


51.9±4.9 


R = </>/K + , <f> 


R = 0.025 ± 0.006 


1.91±0.45 


1.65±0.5 


2.5±0.25 


4.58±0.2 


7.6±1.1 


A 


18.1±1.9 


28±1.5 


41.9±6.1 


43.0±5.3 


44.7±6.0 


44.9±8.9 


A 


0.017±0.005 


0.16±0.03 


0.50±0.04 


0.66±0.1 


2.02±0.45 


3.68±0.55 






1.5±0.13 


2.48±0.19 


2.41±0.39 


3.8±0.260 


4.5±0.20 








0.12±0.06 


0.13±0.04 


0.58 ±0.13 


0.83±0.04 


Q+J2 








0.14±0.07 






K s 












81±4 


V[fm d ] 


3596±331 


4519±261 


1894±409 


1879±183 


2102±53 


3004±1 


T [MeV] 


157.8±0.7 


153.4±1.6 


123.5±3 


129.5±3.4 


136.4±0.1 


136.4±0.1 


\ 


5.23±0.07 


3.49±0.08 


2.82±0.08 


2.42±0.10 


1.94±0.01 


1.74±0.02 




1.657* 


1.41* 


1.36* 


1.30* 


1.22* 


1.16* 


7<7 


0.335±0.006 


0.48±0.05 


1.66±0.10 


1.64±0.04 


1.64±0.01 


1.64±0.001 




0.190±0.009 


0.38±0.05 


1.84±0.32 


1.54±0.15 


1.54±0.05 


1.61±0.02 


A/3 


0.877±0.116 


0.863±0.08 


0.939±0.023 


0.951±0.008 


0.973±0.002 


0.975±0.004 


Mb [MeV] 


783 


576 


384 


344 


271 


227 


Ms [MeV] 


188 


139 


90.4 


80.8 


63.1 


55.9 



tribution. The contamination of A by hyperons decays is 
not material. However, the decays of anti-hyperons con- 
taminates in highly relevant way the yields of p and we do 
not discuss here this effect further, the reader will note the 
relevant yields of p, A and E in table 2. The observable 
yield of p is further contaminated by decays of E + . 

The model yields obtained are shown in table 2. These 
results are prior to any WD contributions. The yields of 
input particles can be compared to the fitted inputs seen 
in table 1. We present also predictions for yields of a num- 
ber of other particles of interest. We do not show the un- 
certainty in these results, which can be considerable: in 
addition to the error propagating through the fit, there is 
systematic error due to the shape of the x 2 /d.o.f. mini- 
mum, see figure 1. 



2.3 RHIC energy range fit 

The RHIC central rapidity particle yields at v / snn = 200 
and 130 GeV are analyzed using nearly the same method 
and principles described in the study of the total parti- 
cle yields. This can be done for the case that the particle 
yields, and hence their source, is subject to (approximate) 
scaling, that means is flat as function of the rapidity dis- 
tribution [55]. The overall normalization of yields then 



contains, instead of the volume V, the volume fraction 
dV/dy associated with the size of the volume at the ra- 
pidity of the source of particles at y. We note that, in the 
local rest-frame, the total yield of particles N± v can be 
written in the equivalent forms: 

[ „r f , dv f , dN 
N 4v = J dVp = J dy—p = J dy—. (3) 



The local rest frame particle density, p 
thus related to the rapidity density by: 

dN _ dV 
dy dy P ' 



dN/dV, is 



(4) 



The SHM fits to particle densities dN /dy thus produce as 
the normalization factor the value dV/dy. The qualitative 
relation between dV/dy and V (rest- frame hadronization 
volume) must include the maximum rapidity range 2y p , 
where y p is the rapidity of the nuclei colliding head on, 



V 



, dV 

k— x 2y 
dy 



pi 



(5) 



where k is a reaction energy dependent constant. The 
study of the total hadron yields at RHIC we present sug- 
gests k ~ 0.4-0.6. 
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Table 2. Output total hadron multiplicity data for AGS (left) and SPS (right). Additional significant digits are presented in 
particle yields for purposes of tests and verification. The statistical parameters generating these multiplicities are the central fit 
values seen in table 1. Hadron yields presented are prior to weak decays and apply to the total multiplicities expected at 
the most central collision bin with the corresponding baryon content b as shown. For SHM parameters, see table 1. 



E [A GeV] 


11.6 


20 


30 


40 


80 


158 


V s nn [GeV] 


4.84 


6.26 


7.61 


8.76 


12.32 


17.27 


Vcm 


1.6 


1.88 


2.08 


2.22 


2.57 


2.91 


centrality 


m.c. 


7% 


7% 


7% 


7% 


5% 


b = B — B 


375.6 


347.9 


349.2 


349.9 


350.3 


362.0 


7T + 


135.2 


181.5 


238.7 


290.0 


424.5 


585.2 




162.1 


218.9 


278.1 


326.0 


461.3 


643.9 


K+ 


17.2 


39.4 


55.2 


56.7 


77.1 


109.7 


IC- 


3.58 


10.4 


15.7 


19.6 


35.1 


54.1 


Rs 


10.7 


25.5 


35.5 


37.9 


55.1 


80.2 




0.46 


1.86 


2.28 


2.57 


4.63 


7.25 


p 


174.6 


161.6 


166.2 


138.8 


138.8 


144.3 


P 


0.021 


0.213 


0.68 


0.76 


2.78 


5.46 


A 


18.2 


29.7 


39.4 


34.9 


42.2 


48.3 


A 


0.016 


0.16 


0.51 


0.63 


2.06 


4.03 




0.47 


1.37 


2.44 


2.43 


3.56 


4.49 


£7 + 


0.0026 


0.027 


0.089 


0.143 


0.42 


0.82 


Q 


0.013 


0.068 


0.14 


0.144 


0.27 


0.38 


n 


0.0008 


0.0086 


0.022 


0.030 


0.083 


0.16 


K°(892) 
A 


5.42 


13.7 


11.03 


12.4 


18.7 


19.1 


38.7 


33.43 


25.02 


26.6 


27.2 


28.2 


A++ 


30.6 


25.62 


22.22 


24.2 


25.9 


26.9 


,4(1520) 


1.36 


2.06 


1.73 


1.96 


2.62 


2.99 


IT (1385) 
~°(1530) 


2.51 


3.99 


4.08 


4.26 


5.24 


5.98 


0.16 


0.44 


0.69 


0.73 


1.14 


1.44 


V 


8.70 


16.7 


19.9 


24.1 


38.0 


55.2 


rf 


0.44 


1.14 


1.10 


1.41 


2.52 


3.76 


P° 


12.0 


19.4 


14.0 


18.4 


32.1 


42.3 


o;(782) 


6.10 


13.0 


10.8 


15.7 


27.0 


38.5 


/o(980) 


0.56 


1.18 


0.83 


1.27 


2.27 


3.26 


(s-s)/(s + s) 





-0.092 


-0.085 


-0.056 


-0.029 


-0.056 



Regarding the data source, and weak decay accep- 
tance, we need to consider case by case the experimental 
results, since the relative importance of hyperon decays in 
the total baryon yields is high. In particular we note: 

— For RHIC-130 dN/dy fit (second column from right in 
the top section of table 3): 

The tt^~ , K^, p and p 5% centrality results are from 
PHENIX [56]. We assume that the Kg decays into pi- 
ons are accepted at 70% level, and Kl at 40% level. 
Mesons (pions and kaons) from hyperon decays are ac- 
cepted at 30% level, while nucleons from hyperon de- 
cay are nearly fully accepted, both 90% and 99% ac- 
ceptances are in essence indistinguishable. 17 ± decays 
are fully accepted. We include in the fit an average 
of the STAR [57] and PHENIX [58] A and A yields 
where we can asses the feed from E and E in view of 
the STAR analysis [59], we accept 99% of Q and Q 
decays into A and A. 

For the E and E weak feed yield corrections are im- 
material. However, we cannot directly use the yields 
as these are presented for the 10% most central reac- 



tions. We fit the weak decay corrected E/A and E/A 
ratios. In order to relate this to the total particle yields, 
we include also E/h~ [h~ = negatives) where we ac- 
cept in STAR h~ the weak decay products according 
to the pattern: Kg decays into pions are accepted at 
90% level, and Kj, at 30% level, pions and kaons from 
hyperon decays are accepted at 50% level, while nu- 
cleons from hyperon decay are accepted, at 99% level. 
The same is assumed in the fit of (12+17)/ h~ also mea- 
sured by STAR [59]. We include in the fit the STAR 
resonance ratios, K°(892)/K" [60] and (j>/K~ [61], in 
both cases we include 50% feed from Q and fl decay 
into kaons, which is immaterial for the result. 
- For RHIC-200: dN/dy fit (last column on right in the 
top section of table 3): 

The tt^~ , K^, p and p 5% centrality results are from 
PHENIX [62] . We assume that the Kg weak decays are 
accepted at 70% level, and at 40% level. Mesons 
(pions and kaons) from Hyperon decays are accepted 
at 30% level, while nucleons from hyperon decay are 
nearly fully accepted, we included this at 90% level 



Jean Letessier, Johann Rafelski: Hadron Production and Phase Changes in Relativistic Heavy Ion Collisions 



7 



Table 3. The input particle data (top) and the resulting statistical parameters, and the chemical potentials derived from 
these, at bottom, for the RHIC energy range. Any of the entries with a * is set as input or is a constraint, e.g., in general A s 
results from the constraint to zero strangeness. * indicates input particle multiplicity derived from interpolating yields between 
different energies, see the RHIC ^snn = 62.4 GeV case. On right, the case of central rapidity yields dN/dy, and on left, the 
total particle yields, in all cases considered for the most central 7% collisions. For N^, we show the participant count. 



V s nn [GeV] 


62.4 


130 


200 


62.4 


130 


200 


E cq [A GeV] 


2075 


9008 


21321 


2075 


9008 


21321 


Ay 


±4.2 


±4.93 


±5.36 


±4.2 


±4.93 


±5.36 




N 47I 5% 


dN/dy\ y=0 5% 


N w 


349±6 


349±6 


349±6 








Q/b 


0.39±0.02 


0.39±0.02 


0.39±0.02 


0.39±0.02 


0.4±0.01 


0.4±0.01 


TV /lX + 








1.02±0.03 


1.0±0.03 


1.0±0.05 


7T+ 


1 1140±90 


1 1450±90 


1677±150 




276±36 


286.4±24.2 


7T~ 






1695±150 




270±36 


001 o i c\r\ a 

281.8±22.8 


trA- 

K + 






293±26 




46.7±8 


48.9±6.3 








243±22 




40.5±7 


45.7±5.2 


0/K 










0.15±0.03 


0.16±0.03 


P 










28.7±4 


18.3±2.6 


P 










20.1±2.8 


I3.5±l.8 


A 








T 17±2 


17.35±0.8 




A 








t 10±l 


12.5±0.8 




S-/hr 










0.0077±0.0016 




S ,5 /A 








t 2.05±0.2 


0.187±0.046 


2.17±0.25 


~ + , ~ + /A 








f 1.3±l 


0.215±0.054 


l.83±0.25 


S /S 










0.853±0.1 




n/h- 










0.0012±0.0005 




(n + ~n)/hr 










0.0021±0.0008 




K°(892)/K" 










0.26±0.08 


0.23±0.05 


V, dV/dy [im] 


4871±394 


6082±384 


8204±351 


932±38 


930±3 


H82±55 


T [MeV] 


140* 


141.9±0.5 


142.4±0.01 


142.2±0.01 


143.8±0.1 


I4l.5±0.l 


\ 


1.35±0.02 


1.25±0.01 


1.20±0.01 


1.15±0.02 


1.076±0.001 


l.062±0.00l 


Xs 


1.104* 


1.074* 


1.069* 


1.054* 


1.025* 


1.024* 


7<7 


1.62* 


1.62* 


1.62* 


1.62* 


1.59±0.001 


l.56±0.0l 


Is 


2.18±0.2 


2.20* 


2.00±0.29 


2.13±0.14 


2.22±0.01 


2.00±0.02 


A/3 


0.933±0.001 


0.979±0.001 


0.988±0.002 


0.986±0.002 


0.997±0.001 


0.997±0.00l 


/is [MeV] 


126 


94.8 


79 


61.2 


31.5 


25.7 


/is [MeV] 


27.7 


21.4 


16.5 


13.6 


7.0 


5.2 



in the reported fit. S ± decays are fully accepted. We 
take STAR resonance ratios, K°(892)/K- [63,64] and 
<fi/K~ [65], in both cases we include 50% feed from 
Q and fl decay into kaons. The method to study the 
yields of stable hadrons along with resonances follows 
the work on the impact parameter dependence at ^/snn 
= 200 GeV [3]. 

We did not use yields of A and A since without direct 
measurement of S and S it is hard to judge the weak 
decay contamination in the data. Furthermore, we pre- 
ferred to study the relative yields p/ir + ,p/ir~ . In the 
fit presented, we assumed that the pion feed from WD 
of hyperons is at 80% level. The other WD character- 
istics are as discussed just above. This slight change in 
data input and also the slight modification of the pat- 
tern of weak decay acceptance has, in comparison to 
Ref . [3] , yielded a increase of the volume factor dVj dN 
by 1.2 s.d., while other variations are within 0.5 s.d.. 



We can expect, in near future, particle multiplicity re- 
sults from RHIC obtained at ^/snn=62.4 GeV. We in- 
terpolate the central rapidity yields of strange hyperons 
A, A, E and S, presented in [66], to this energy. With 
these 4 inputs, two constraints, setting the 7 9 = 1.62, 
T = 140, we find a good description of the interpolated 
data but with a few degrees of freedom. We have four 
interpolated 'data' points, two constraints — strangen- 
ess conservation and Q/b, thus 6 data points which are 
fitted using four flexible parameters, T,dV/dN, \ q , and 
A/ 3 . This set of parameters, then, yields our prediction of 
central rapidity particle multiplicities, seen in table 4, for 
yii5"=62.4 GeV. 

We make an effort to understand also the recently 
finalized total multiplicities of and n ± [67] at 
yl snn = 200 GeV. Additional qualitative constraint is de- 
rived from total charge particle multiplicities [68] , however 
this result is not used directly in the fit. With the three 
constraints, four BRAHMS particle yields, we have 7 data 
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Fig. 1. xVdof (top) and the associated significance level P[%] 
(bottom) as function of 7 9 , the light quark phase space oc- 
cupancy, for the AGS/SPS energy range and for the (central 
rapidity) RHIC results. 



points, and also 7 SHM parameters. To be able to make 
a fit with at least one degree of freedom it is necessary to 
make some 'natural' hypothesis. For this reason, we do not 
discuss the fit quality of yields at RHIC but we dis- 
cuss the expected total particle yields, which we regard to 
be an experimentally motivated hadron yield prediction. 

We choose to consider 7 9 = 1.62 ~ e m 7r /2T ^ we 
find systematically at the RHIC energy scale. Our 'fit' to 
Ni„ data at ^/snn = 200 GeV works, but it must not be 
seen as a full fit, rather a consistency test of SHM. This 
is allowing a prediction to be made of other AT 47r we show 
in table 4. This consideration is also yielding a rapidity- 



averaged value of T and of the 5 chemical parameters, as 
well as an estimate of the proper size V of the hadronizing 
fireball. The value of hb, which varies as function of rapid- 
ity, following the highly variable baryon distribution [69], 
is found at a median value, seen at the bottom of table 3, 
on left for the N± v fits. 

We extend the consideration of the N± v yields to the 
lower energies, v /snn = 62.4 and 130 GeV. This can be 
done assuming that there is no change in physics between 
top SPS energy and RHIC 200 GeV run. Thus, the success 
of our particle yield prediction would be a confirmation of 
this hypothesis. Our procedure can be seen in detail on 
the left of table 3. We fix the hadronization temperature 
at T = 140 MeV, choose the value j q = 1.62 ~ e m ^ 2T , 
and interpolate the values of 7 S . We do find the required 
values of X q , A/3 and V needed to assure the total baryon 
yield, fraction of charge Q/b and one particle yield, which 
we choose to be the interpolated total ir + . We use the 
observation that the ir + yield from Brahms connects, in 
a logarithmic plot, in a nearly perfect straight line with 
the SPS energy domain. This produces the 7r + interpo- 
lated values we introduced in table 3. The SHM succeeds 
perfectly and allows us to offer predictions for the total 
particle yields presented in table 4. 

We present, in detail, the resulting particle multiplic- 
ities in table 4 for RHIC. On left, we show the expected 
total yields and on right the central rapidity yields. We re- 
call that, among total yields, only at 200 GeV a significant 
experimental input was available, thus the 62.4 and 130 
GeV total yield results are an educated guess satisfying all 
constraints and criteria of the SHM model. Similarly, the 
central rapidity region yields for 62.4 GeV is a prediction 
based on interpolated yields, with inputs seen in table 3. 
All results, presented in table 4, are obtained prior to WD. 



2.4 Energy dependent particle yields 

We consider, more systematically, the energy dependence 
of particle yields and ratios. Of particular interest is the 
ratio K + /-7r + which shows the previously unexplained horn 
structure. We compare the experimental and theoretical 
behavior in figure 2. The 47r results are blue filled squares. 
The central rapidity RHIC results (on right in red) are 
shown as open squares, while the predicted total yield 
ratio for ^/snn = 62.4 is given as an open circle. 

We recall that the abrupt increase in the value of j q 
occurs where the rise in K + /7r + reverses, turning into 
a sudden decrease with reaction energy. The solid line, 
shows our chemical non-equilibrium fit which reproduces 
the horn structure well. The predicted total yield ratios 
for v / snn = 62.4 and 130 GeV (edges in solid line) arise 
from the interpolation of yields and/or continuity in value 
of statistical parameters such as 7 g between the top SPS 
and the top RHIC energy, see above, subsection 2.3. 

The dotted line, in figure 2, presents best fit results 
obtained within the chemical equilibrium model, i.e., with 
7« = 7g = 1) using the same computer program (SHARE), 
and the same data set. We see that the chemical equilib- 
rium SHM cannot explain the horn in the K + /n + ratio. 
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Table 4. Output hadron multiplicity data for the RHIC energy range. See text for the meaning of predictions of yields 
at 62.4 and 130 GeV and of dN/dy at 62.4 GeV. The input statistical parameters are seen in table 3. b = B — B = Nw for 47r 
results and b — d(B — B) /dN for results at central rapidity. Additional significant digits are presented for purposes of tests and 
verification. All yields are without the weak decay contributions. 



\/ -snn [GeV] 


62.4 


130 


200 


62.4 


130 


200 


£ eq [GeV] 


2075 


9008 


21321 


2075 


9008 


21321 


Ay 


±4.2 


±4.93 


±5.36 


±4.2 


±4.93 


±5.36 




JV 4t 5% 


dN/dy\ y=0 


5% 


b 


350.2 


350.2 


349.6 


33.48 


18.50 


14.8 


7T + 


899 


1201 


1543 


183.8 


230.3 


239.8 


7T~ 


927 


1229 


1573 


186.7 


231.9 


241.0 


K+ 


230.9 


302.5 


291.9 


43.7 


47.9 


47.1 


KT 


168.5 


238.4 


242.3 


37,6 


44.2 


44.2 


K S 


193.8 


261.0 


259.9 


39.4 


44.4 


44.2 


<t> 


27.3 


34.6 


28.9 


5.74 


6.86 


6.18 


p 


140.0 


157.6 


192.0 


19.34 


17.09 


16.34 


P- 


24.1 


42.9 


66.1 


8.37 


11.11 


11.44 


A 


81.1 


97.4 


89.9 


12.3 


12.04 


10.7 


A 


20.2 


35.1 


38.3 


6.36 


8.60 


8.02 




12.9 


16.4 


11.6 


2.14 


2.30 


1.91 




4.6 


7.79 


6.13 


1.32 


1.80 


1.53 


Q 


1.94 


2.68 


1.45 


0.36 


0.44 


0.33 


n 


1.04 


1.74 


0.98 


0.27 


0.38 


0.29 


K°(892) 
A 


48.7 


67.4 


68.1 


10.2 


11.5 


11.2 


27.6 


31.1 


38.1 


3.78 


3.32 


3.15 


A++ 


26.2 


29.9 


36.9 


3.69 


3.30 


3.13 


71(1520) 


4.43 


6.4 


6.0 


0.81 


0.81 


0.70 


r+(1385) 


9.80 


11.91 


11.19 


1.52 


1.50 


1.33 


=■0(1530) 


4.20 


5.46 


3.88 


0.71 


0.78 


0.64 


V 


131.9 


179.5 


192.3 


27.2 


30.5 


30.6 


v' 


10.8 


15.2 


14.64 


2.30 


2.64 


2.51 


P° 


85.8 


117 


157 


18.1 


19.5 


20.3 


^(782) 


75.9 


104 


142.8 


16.2 


17.4 


18.3 


/o(980) 


6.51 


9.03 


12.96 


1.40 


2.02 


1.58 


(a - a) /{a + a) 





















The dashed line corresponds to the result obtained fixing 
7 9 = 1 but allowing 7 S to assume a best value. We see that, 
without 7 9 > 1, it is difficult if not impossible to obtain 
the large reduction of K + /7r + ratio with increasing energy. 
These findings are in line with prior attempts to explain 
the horn-structure, see, e.g., figure 4 in [54]. We note that 
our semi-equilibrium results follow better the trend set by 
the experimental data, which is a consequence of the re- 
laxation of exact strangeness conservation requirement. It 
appears that the full chemical non-equilibrium statistical 
hadronization model is required in order to obtain sat- 
isfactory understanding of the energy dependence of the 
K+/7T+ ratio. 

A graphic comparison of the experimental input, and 
theoretical output particle yields as function of energy for 
several other particles is seen in figure 3. We show K~ /n~ 
with A/ir+, ~-/n - with 5Q-/w-, E + /A with 2l2/E + 
and at bottom (j>/V tt + tt~ . We are showing the total SHM 
yield ratios at AGS/SPS as well as at RHIC, connected 
by a solid line (thick and/or thin). The central rapidity 



yields at RHIC are also presented for comparison by the 
dashed lines, (in red and violet) . 

The SHM allowing for chemical non-equilibrium repro- 
duces all salient features of the experimental particle yield 
data well as function of energy, including the NA49 results 
that otherwise could not be described in equilibrium and 
semi-equilibrium approach [22], e.g., the already discussed 
K + /n + ratio shown in figure 2. In addition, In figure 3, 
we note in the top panel the shift of s quark population 
from its dominant baryon component (see A/tt + ) at low 
reaction energy, to meson carriers (see K~ /ir~). 

Of particular importance, in the study of quark-gluon 
plasma formation, is strange antibaryon enhancement. It 
is one of important signatures of deconfinement [70] . These 
particles are hard to make in conventional environment, 
and also are highly sensitive probes of the medium from 
which they emerge. There is still only fragmentary data 
available for antibaryon ratios of interest, EE + / A and 2fi/E + , 
shown in figure 3, in the third panel from the top. 

In figure 4, we show as function of energy the individ- 
ual yields of A, S + and Q + , along with the experimental 
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0.2 



0.15 




0.2 



Fig. 2. K + /n + total yields (filled squares, blue) and cen- 
tral rapidity density (open squares, red) as function of snn • 
The solid lines show chemical non-equilibrium model fit. The 
chemical equilibrium fit result is shown by the dotted line. 
The dashed line arises finding best y 3 for ^y q — 1. See text 
about the total yield results at ^/snn = 62.4 and 130 GeV 
(unmarked edges in lines) and about the central rapidity yield 
at Snn = 62.4 (open circle). 



data used as input. All three antihyperon production rates 
are predicted to rise at nearly the same rate as function 
of a/snn — indeed the strange antibaryon ratios, we have 
shown in figure 3, are as good as flat compared to the great 
variability of the absolute yields, which are increasing very 
rapidly, as is seen in figure 4. The strange antibaryon pro- 
duction enhancement phenomenon has been considered in 
terms of a comparison of yields of antibaryons to a baseline 
yield obtained scaling the pp or pBe yields. This was done 
as function of impact parameter, and reaction energy [71]. 

We note that the SHM non-equilibrium approach undcr- 
predicts by 1.5 s.d. the yields of Q and f2, for both SPS- 
NA49 results available at 40 and 158 A GeV. We did not 
include 158 A GeV results in our input data set since non- 
SHM source, such as chiral condensate [23], can generate 
such an excess, 

Ratios of strange antibaryons are a sensitive probes 
of the medium. Once the deconfined phase is formed, the 
ratios of yields of strange antibaryons should not change 
drastically. Thus, our finding as function of energy in essence 
of a flat ratio, with minor fluctuations originating in the 
other experimental data and amplified by the sensitivity 
of these particles, suggest that the same form of (decon- 
fined) matter is present at SPS and RHIC, except perhaps 
for the lowest SPS reaction energy. 

Another important point is that these antihyperon ra- 
tios are relatively large, and hard to understand except in 
terms of the quark coalescence picture. It would be very 
interesting to confirm experimentally that, at AGS energy 

scale, 2{]/S + is indeed as large as predicted in figure 3. 
This would establish coalescence quark chemistry in this 
low energy environment. Further, this maybe taken as an 
indication that the transition we observer at ,AS N in- 




10 100 
Vs^J [GeV] 



Fig. 3. Comparison of experimental and theoretical ratios of 
particle yields as function of reaction energy ^/snn — theoret- 
ical SHM total Nttn results are connected by solid lines, with 
thick/thin lines corresponding to the different particle ratios. 
The (corresponding) dashed lines connect the central rapidity 
dN/dy results at RHIC. Experimental data used as fit input is 
shown with its error bar. 



volves two deconfined phases of different structure. We 
will further discuss this in section 4. 



2.5 Yields of pentaquark hadrons 

There are now more than 600 papers with the title word 
'pentaquar', however, on balance the evidence for the ex- 
otic hadrons <9 + (1540) with the quark content [uudds], 
and typical decays pK° ,nK* + , and S (1862) [ssddu], 
with typical decays 5~it~ , , is not convincing. The 

a (1862) may have been observed by NA49 in pp in- 
teractions at top SPS energy [72]. The <9 + (1540) remains 
uncertain: several high statistics confirmation experiments 
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[GeV] 



Fig. 4. Yields of strange antibaryons as function from 
top to bottom A, E + and fi + . The solid lines connect the re- 
sults of SHM N±tt fit to particle data. The AGS/SPS energy 
range yields (in blue) on left and RHIC N± n on right. Also 
on right (in red, connected with dashed lines) are the central 
dN/dy yields. The yields at ^/snn = 62.4 used in our study 
are result of interpolation of RHIC and SPS results. 



failed to find this state. Arguments were presented why 
these states were incorrectly identified [73] . We have there- 
fore not included these and other related exotic quark 
states from the hadron resonance list in SHARE, when 
performing the fits here presented. 

On the other hand, theoretical arguments for the ex- 
istence of pentaquark states have not been refuted. We 
thus present in table 5 predictions for the production rates 
of + and E~~ . These yields are highly sensitive to the 
hadronization conditions [74] , and were obtained using the 
parameters of the fits here presented. 

The expected statistical hadronization yield of <9 + rises 
rapidly, by an order of magnitude, between 11 and 30 
A GeV reaction energies and remains practically constant 
thereafter. The expected (9 + (1540) yield in fact exceeds 
the SHM predicted yield of ^1(1520) in the threshold en- 
ergy domain by a factor 2-4, and comparing to the ob- 
served /1(1520) yield at 158 A GeV by more than a factor 
4. Furthermore, at 30 A GeV the background multiplic- 
ity is relatively small, while the rapidity range is also re- 
stricted compared to the top SPS energy, which should 
help finding the pentaquark, if it exists in the range of 
energies characterizing the horn in the K + /ir + . 



3 Fireball properties at breakup 

3.1 Energy dependence of model parameters 

The statistical parameters of the SHM are shown, as func- 
tion of y/sNN, in figure 5, for the entire energy domain. 
From top to bottom, we see the chemical freeze-out tem- 
perature T, the statistical occupancy parameters 7, and 



170 
,— , 160 
% 150 
<-> 140 
^ 130 

120 

1.5 



a* 



1 

0.5 






[GeV] 



Fig. 5. Statistical parameter results for (blue online, 

square). From top to bottom: T, 7,, 7 9 /7s and (j,b , [MeV] , 
as function of ^/snn- The lines guide the eye. Same for dN/dy 
at RHIC (red online triangles). 



7 S /7 9 and the chemical potentials /xb and fig. The error 
bars comprise the propagation of the experimental yield 
errors, as well as any uncertainty due to the shape of the 
X 2 /d.o.f. minimum, seen in figure 1. The (red) triangle 
results are for the RHIC dN/dy case, while (blue) squares 
are for the data throughout the energy domain and 
include the estimates we made for the RHIC energy range 
(at yj snn = 62.4 and 130 GeV, we do not show for these 
fits an error bar, as these results are solely our estimate). 

The only significant difference between RHIC dN/dy 
and N^n results is noted for the chemical potentials /^b 
and (is and shown in the bottom panel (note logarithmic 
scale). The baryochemical potential /j,b drops relatively 
smoothly as the reaction energy is increased. The vertical 
line indicates the observed sudden change in the structure 
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Table 5. Predicted yields of ©+(1540) and 5~ 
shown in tables 1 and 3. 



"(1820) pentaquarks for AGS, SPS and RHIC, obtained with SHM parameters 



E [A GeV] 


11.6 


20 


30 


40 


80 


158 


2075 


9008 


21321 


[GeV] 


4.84 


6.26 


7.61 


8.76 


12.32 


17.27 


62.4 


130 


200 


centrality 


m.c. 


7% 


7% 


7% 


7% 


5% 


5% 


5% 


5% 


(9+ (1540) 


0.66 


1.14 


6.90 


7.15 


6.52 


6.70 


7.23 


7.92 


7.19 


=^-(1820) 


0.0022 


0.010 


0.098 


0.11 


0.18 


0.24 


0.71 


0.89 


0.78 



of the fireball. This is seen in all statistical variables, but 
most clearly in j q . 

It is important to recall that we present 7, evaluated 
using hadronic multiplicities. If these arise from breakup 
of a quark fireball, the quark-side occupancy parameters 
could be considerably different. The hadron-side phase 
space size is, in general, different from the quark-side phase 
space, since the particle degeneracies, and masses, are 
quite different. In the study of the breakup of the quark 
fireball into hadrons, we can compute the resultant hadron 
phase space occupancy for two extreme limits. 

First, consider a fast transformation (sudden breakup) 
of the quark phase. This occurs nearly at fixed volume. To 
accommodate the difference in the momentum part of the 
phase space, the chemical occupancy non-equilibrium pa- 
rameters ji undergo an abrupt change. We note that it is 
of no importance if there was or not a phase transition be- 
tween the phases, what matters is that there was no time 
to rccquilibrate chemically the quark yields. In the oppo- 
site limit of a very slow transformation of phases, there is 
available a prolonged period in time in which the volume 
of the system can change to accommodate the appropriate 
number of particles in chemical equilibrium corresponding 
to the maximum entropy content. 

To determine the change in 7; in sudden hadroniza- 
tion, one needs to compare in detail the phase space of 
quark phase with that of hadron gas. In order to make 
this comparison one must consider the energy and en- 
tropy content of the QGP phase. For /iB = 0, as well 
as small chemical potentials hb/T < 1, lattice evaluation 
of the deconfined phase properties are available [75-78]. 
It is thus possible to model quantitatively the properties 
of the deconfined phase, and to compare these with the 
results of the SHM [79]. The remarkable result is that 
near to T = 140 MeV the sudden transition requires the 
value 7 9 — 1.6 on the hadron side, if the quark phase was 
chemically equilibrated. 

Thus, the large values of 7^° and 7^ G , seen in figure 5 
at large ^/snn, where hb/T < 1 arc consistent with sud- 
den breakup of chemically equilibrated primordial QGP 
phase. Other dynamic effects, in particular fast expan- 
sion, in general also favor an over-saturated phase space 
with 7j > 1. 

As seen in figure 5, j s /j q rises at first rapidly, as ex- 
pected if strangeness production is delayed by a greater 
threshold mass and has to catch up with the light hadron 
production. "f s /"f q decreases beyond the edge of energy 
threshold, as can be expected due to the conversion of the 
quark to hadron occupancy discussed above. The rise re- 



sumes and continues for all energies above 80 AGeV (note 
that, at RHIC only, results showing an error can be con- 
sidered to arise from a fit). 



3.2 Physical properties 

We now turn our attention to the physical properties of the 
hadronizing fireball obtained summing individual prop- 
erties of hadronic particles produced. One can view the 
consideration of the physical properties of the fireball at 
breakup as another way to present the SHM parameters. 
For example, the net baryon density, pb = (B — B)/V, is 
most directly related to the baryochemical potential /xb, 
the thermal energy density E t h/V is related to T etc. 

We present the physical properties, i.e., pressure P, 
energy density e, entropy density S/V, net baryon density 
Pb = (B — B)/V and the yield of strangeness s, in table 
6 for the AGS/SPS energy range considered. Note that 
s contains hidden strangeness from 77, <j> and rj'. At the 
bottom of table 6, we show the dimensionless ratios of 
extensive variables P/e, and E t h/TS. These two ratios 
are very smooth as function of energy, and lack any large 
fluctuations that could be associated with fit error. These 
ratios are characteristic for the conditions of the fireball 
at the point of hadronization. 

The results presented can be used to constrain dy- 
namical models describing the evolution of the QGP fire- 
ball in time towards hadronization/particle freeze-out. We 
present the energy range at RHIC on left in table 7. We 
recall that the 62.4 GeV and the 130 GeV 4tt results, as 
well as in part the 200 47T results, are result of consider- 
ations which do not involve experimental measured par- 
ticle yields. Thus, the 47r results are to be seen as SHM 
sophisticated prediction. On right, in table 7, we present 
the results for central rapidity densities. Here, only the 
62.4 GeV case is a prediction, the other results are direct 
consequence of the data interpretation in terms of SHM. 

The fit uncertainty in the quantities presented in ta- 
bles 6 and 7 is difficult to evaluate in detail. The indi- 
vidual physical properties require powers and exponents 
of statistical parameters, and thus, at first sight, we ex- 
pect that the fractional errors are increased, as compared 
to those prevailing among statistical parameters in table 
1. However, the dominant contributions to each physical 
property is often directly derived from the individual ob- 
served particle yields. Therefore, a large compensation of 
errors originating in the fitted statistical parameter errors 
must occur. 
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Table 6. The physical properties: Pressure P, energy density e = E t h/V, entropy density S/V, strangeness density s/V for 
AGS and CERN energy range at, (top line) projectile energy E [GeV]. Bottom: dimensionless ratios of properties at fireball 
breakup, P/e and E th /TS . 



E[AGeV] 


11.6 


20 


30 


40 


80 


158 


V s nn [GeV] 


4.84 


6.26 


7.61 


8.76 


12.32 


17.27 


P[MeV/fm 3 ] 


21.9 


21.3 


58.4 


68.0 


82.3 


76.9 


e[MeV/fm 3 ] 


190.1 


166.3 


429.7 


480.2 


549.9 


491.8 


5/\/[l/fm 3 ] 


1.25 


1.21 


2.74 


3.07 


3.54 


3.26 


100s/V[l/fm 3 ] 
p 6 [l/fm 3 ] 


0.988 


1.52 


5.32 


5.85 


7.65 


7.24 


0.104 


0.0753 


0.184 


0.186 


0.167 


0.121 


P/e 


0.115 


0.128 


0.136 


0.142 


0.150 


0.156 


E th /TS 


0.96 


0.92 


1.27 


1.20 


1.14 


1.11 



Table 7. The physical properties for RHIC energy range, see table 6 for details. For the central rapidity case, we show the 
rapidity densities: energy rapidity density de = dE t h/dV, entropy rapidity density dS/dV, strangeness rapidity density ds/dV 
and net baryon rapidity density db/dV. All 62.4 GeV results, and the 130 GeV 4ir result are, as discussed in text, result of 
assumptions, and/or interpolations of yields and/or parameters, and hence are a prediction. 



V s nn [GeV] 


62.4 


130 


200 


62.4 


130 


200 






dN/dy\ y= 





P[MeV/fm 3 ] 


82.4 


87.8 


80.0 


80.5 


91.4 


94.5 


dE th /dV[MeV/fm :i ] 


516.6 


548.4 


478.9 


532.5 


604.4 


479.4 


dS/dV[l/im 3 ] 


3.62 


3.73 


3.32 


3.64 


4.03 


3.32 


100ds/dV[l/fm 3 ] 
W0db/dV[l/fm 3 ] 


11.5 


12.4 


9.2 


12.0 


13.7 


10.4 


7.19 


5.76 


4.26 


3.59 


1.99 


1.26 


PdV/dE th 
dE th /TdS 


0.159 
1.02 


0.160 
1.04 


0.167 
1.01 


0.151 
1.03 


0.151 
1.04 


0.197 
1.02 



For example, most of the pressure at breakup is due to 
the most mobile, lightest particle, the pion. These yields 
are known to better than 10%, and thus, the pressure must 
be known to greater precision since there are further con- 
straints from consistency of this yield with the yield of 
other particles. This explains why the results when pre- 
sented graphically (see figure 6) are at 5% level smooth 
functions of ^/snn, with fluctuations apparently at worse 
similar to those we see in the individual statistical SHM 
parameters. In future, one could hope to fit to the exper- 
imental data directly the physical properties, bypassing 
the statistical parameters. This can be done, in principle, 
considering the mathematical properties of these expres- 
sions. However, such study transcends considerably the 
scope of this paper, and it is indeed motivated by results 
we present for the first time here. 

On the left hand side, in figure 6, we see from top to 
bottom the baryon density, the thermal energy density e 
and the entropy density a. On the right hand side, from 
top to bottom, we show the pressure P, and the dimen- 
sionless ratio of pressure to thermal energy density P/e, 
and Eth/TS = e/Ta. The triangles (red) correspond to 
the properties of the fireball at central rapidity at RHIC 
energy scale. We note a significant difference between the 
total fireball averages (squares) and the central rapidity 
result (triangles) only in the net baryon number density. 



As the reaction energy passes the threshold, 6.26 GeV 
< ^/sJra < 7.61 GeV, the hadronizing fireball becomes 
much denser: the entropy density jumps by factor 4, and 
the energy and baryon number density by a factor 2-3. 
The hadron pressure increases from P = 25 MeV/fm 3 
initially by factor 2, and ultimately more than factor 3. 
There is a more gradual increase of P/e = 0.115 at low 
reaction energy to 0.165 at the top available energy. Also 
Et h /TS falls gradually from 0.9 down to 0.78 for the high 
density fireball. 

The rather rapid change in the individual properties: 
entropy, energy, pressure is seen, in figure 6, to be largely 
compensatory, resulting in a smooth change in P/e, and 
similarly Eth/TS. Even though there is a small residual 
variation reminding us of the sudden changes in the three 
factors in the ratio E/TS, this quantity is extraordinarily 
smooth. Moreover, we see the same value for the central 
rapidity as we find for the average over the entire fireball. 
Thus E/TS ~ 0.78 could be a universal hadronization 
constraint. 

For AGS 11 GeV and SPS 20 GeV results, the value 
E/TS is greater, reaching to E/TS ~ 0.9. This requires 
P/e to be smaller, as is the case when the effective quark 
mass increases. A simple structure model employing a 
thermal quark mass, m q ~ aT, was considered in Ref. [3]. 
E/TS ~ 0.78 corresponds to the value of a ~ 2 (usual for 
thermal QCD) and found in the limit of large A and here 
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Fig. 6. From top to bottom, on left hand side the baryon density p B _-g[fm 3 ], energy density e[MeV/fm 3 ], and entropy density 
a-[fm~ 3 ], as function of ^/snn, on right hand side pressure P[ MeV/fm 3 ], P/e, and E th /TS = e/Tcr. Squares (blue) average over 
the entire fireball at hadron freeze-out, triangles (red) for the central rapidity region of the fireball. 



large ^/snn- a has to rise to a > 4 in order to explain the 
rise in E/TS. This points to a phase of 'heavy' quarks be- 
ing at the origin of the increase of E/TS with decreasing 
^/snn- Such a heavy quark phase can be possibly a 'valon' 
quark phase and we pursue this further in section 4.4. 

We believe that any structure model of the phase trans- 
formation, and/or the two phase structure will need to ad- 
dress E/TS, and/or P/e freeze-out condition results quan- 
titatively. These two ratios, E t h/TS and P/e, are related. 
Restating the 1st law of thermodynamics: 

+ = ! + _!_! i, k = —. 6 

IS a e 

For each hadron fraction with density the total fugacity 
is 

li;v/. (7) 

3 

where all valance quarks and antiquarks of each hadron 
fraction contribute in the product, see section 2 in [14]. In 
the limit of chemical equilibrium: 



a To 



(8) 



Thus, in this limit at the RHIC energy range, we expect 
that Eth/TS — > l/(l+P/e). However, the results in figure 



6 show, the chemical non-equilibrium effects contribute 
considerably. 

It is interesting to note that the same behavior of the 
physical properties of the fireball has also been obtained as 
function of the volume in the study of impact parameter 
dependence, see figure 4 in [3]. In fact, the results we de- 
rived show an unexpected universality of the hadronizing 
fireball, which depends solely on the question if it occurs 
'below' or 'above' the threshold in energy and volume size; 
the volume threshold corresponds to critical participant 
number 13.4 < 4 cr < 25.7. 

At these values of A and the associated baryon content 
at central rapidity [3], the grand canonical description of 
particle yields is still justified [80,81], also for strangeness. 
However, the fitted reaction volume (not further used in 
the present work) may be revised within the canonical 
approach by 10 - 20% for the most peripheral collisions 
studied in [3]. 



3.3 Strangeness and Entropy yield 

The yield of strangeness produced, should the deconfined 
QGP fireball be formed, is sensitive to the initial con- 
ditions, especially to the temperature achieved. The stan- 
dard results for strangeness relaxation time corresponds to 



Jean Letessier, Johann Rafelski: Hadron Production and Phase Changes in Relativistic Heavy Ion Collisions 



15 



t s (T = 300 MeV) ~ 2fm/c [82]. When this result is used 
in model calculations addressing RHIC [83], one finds, as- 
suming gluon thermal and chemical equilibrium, that the 
thermal strangeness production in the early stage suffices 
to saturate the QGP fireball phase space at hadronization. 

Even so, there is considerable uncertainty how short 
the time required to relax strangeness flavor is, as the re- 
laxation time lengthens with the square of the glue phase 
space under-occupancy, t s oc 1/7q- Much of the uncer- 
tainty about the gluon chemical conditions prevailing in 
the initial thermal phase can be eliminated by consider- 
ing the ratio of the number of strange quark pairs to the 
entropy s/S. In the QGP, the dominant entropy produc- 
tion occurs during the initial glue thermalization 70 — ► 1, 
and the thermal strangeness production occurs in parallel 
and/or just a short time later. Moreover, both strangen- 
ess s and entropy 5* are nearly conserved in hadronization, 
and thus, the final state yield value for the ratio s/S is di- 
rectly related to the kinetic processes in the fireball at 
r ~ 1-3 fm/c. A thorough discussion of the observable 
s/S is presented in [79], and detailed evaluation within a 
dynamical model of s/S was obtained [84]. The following 
is a motivating introduction to these developments. 

We first estimate the magnitude of s/S in the QGP 
phase considering, in the hot early stage of the reaction, an 
equilibrated non-interacting QGP phase with perturbative 
properties: 

£ = P± _ 7s QGP (3A 2 )T 3 (m s /T) 2 K 2 (m,/T) 
S ~ a (32tt 2 /45)T 3 + n f [(7tt 2 /T5)T 3 + /i 2 T] ' 

l + 0.054(lnA 9 ) 2 ' v ' 

Here, we used for the number of flavors rif = 2.5 and 
m s /T = 1. We see that the result is a slowly changing 
function of X q , for large \ ~ 4 we find at lowest SPS en- 
ergies, the value of s/S is reduced by 10%. Considering the 
slow dependence on x — m s /T ~ 1 of W(x) — x 1 Ki{x), 
there is further dependence on the temperature T. 

The rise with reaction energy toward the limiting value, 
s/S — 0.03 for large ^/snn> is driven by the decrease 
in X q — > 1 and, importantly, by an increase in chem- 
ical strangeness equilibration with the QGP occupancy 
7? GP — * 1- The dependence on the degree of chemical 
equilibration which dominates the functional behavior with 
V s nn is: 

£_ 0.03 7 Q GP 

S ~ 0.387G+ 0.12 7s QGP + 0.5 79 QGP + 0.054 7 , QGP (ln \ q f ' 

(10) 

Eq. (10) predicts a smooth increase in s/S toward its 
maximum value which by counting the degrees of freedom 
appears to be s/S — > 0.03, while the QGP source of parti- 
cles approaches chemical equilibrium with increasing col- 
lision energy and/or increasing volume. It is important to 
keep in mind that the ratio s/S is established early on in 
the reaction, and the above relations and associated chem- 
ical conditions we considered apply to the early hot phase 
of the fireball. Yet, strangeness and entropy, once created, 



cannot disappear as the more complex low temperature 
domain is developing. Specifically near to hadron freeze- 
out, the perturbative QGP picture used above does not 
apply. Gluons are likely to freeze faster than quarks and 
both are subject to much more complex non-perturbative 
behavior. But the value s/S is preserved across this non- 
perturbative domain. 

In tables 8 and 9, we present, in top portion, the strang- 
eness production as function of reaction energy at AGS, 
SPS and RHIC, respectively. We give the baryon content 
and the total strangeness content of the fireball derived 
from the SHM fit to particle yield. Below, we see the above 
discussed strangeness per entropy s/S ratio, and strang- 
eness per net baryon number s /b ratio. We present the in- 
creasing specific strangeness per baryon and per entropy 
yields in figure 7, two top panels. The remarkable result 
we find is that the specific per entropy yield of strangen- 
ess converges for top RHIC energy and central rapidity 
toward the QGP result obtained counting the degrees of 
freedom, see Eq. (9). The somewhat smaller values for the 
47r case are consistent with the average being made over 
the fragmentation region. This effect is greater in the ratio 
s/b as we have to count all participant baryons. 

The middle section, in table 1, shows the center of 
momentum energy cost ^/snn/(2s/6) to make one strang- 
eness pair. The micro canonical input variables, s/V and 
b/V, for this entry vary significantly along with, and as 
function of ^/%u. Yet, we see that the result obtained 
varies smoothly, at first it diminishes finding a minimum 
at around E — 40 A GeV and it rises slowly thereafter. It 
is clearly more energy expansive to make strangeness at 
AGS, nearly by factor 2. A minimum in energy cost to 
make strangeness is near to 30 ^4GeV beam energy, at the 
peak of K+/7T+ horn. 

The increase in cost of making strangeness can be at- 
tributed to the decreasing energy fraction stopped in the 
reaction. The energy stopping can be estimated by eval- 
uating the per baryon thermal energy content E t h/b and 
obtaining from this the fraction of the initial energy con- 
verted into thermal energy in the final state, (2E t h/b)/ 
\J snn , which fractions steadily drops from 75% at AGS 
to 48% at top SPS energy. 

In terms of thermal energy, the cost of making strang- 
eness pair is given in the last line of table 1. After an 
initial very rapid drop from AGS cost at 20 GeV to 8 
GeV near to the top of the horn, there follows a very 
slow and gradual decrease. We show this result graphi- 
cally in the bottom panel in figure 7. This behavior clearly 
shows a rapid but smooth change-over in the underlying 
mechanism of strangeness production with increasing re- 
action energy, between 11.6 and 30 AGeV. Once the new 
mechanism is fully operational, we have essentially a flat, 
slowly decreasing energy cost per strangeness. The drop 
we observe above 30 AGcV can be thought to originate 
in transfer of thermal energy to the kinetic energy of col- 
lective expansion which we do not record in our analysis, 
and thus, it is conceivable that the cost in actual energy 
remains constant above ,/s™. 
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Table 8. AGS and CERN energy range (see top lines for projectile energy E [GeV] and v / 3nn): Strangeness yield s (= s), 
strangeness per entropy s/S, strangeness per baryon s/b, the energy cost to make strangeness pair ^/snn/(2s/&), thermal energy 
per baryon at hadronization E th /b, fraction of initial collision energy in thermal degrees of freedom, (2_E t h/6)/V s NN- 



E [A GeV] 


11.6 


20 


30 


40 


80 


158 


V s nn [GeV] 


4.84 


6.26 


7.61 


8.76 


12.32 


17.27 


b 


375.5 


347.9 


349.2 


349.9 


350.3 


362.0 


s 


35.5 


70.3 


100.8 


110 


161 


218 


WQs/S 


0.788 


1.26 


1.94 


1.90 


2.16 


2.22 


s/b 


0.095 


0.202 


0.289 


0.314 


0.459 


0.60 


y/s^/(2S/b) [GeV] 


25,5 


15.5 


13.1 


13.9 


13.4 


14.4 


E th /b [GeV] 


1.82 


2.26 


2.33 


2.58 


3.30 


4.08 


(2£ t h/b)/y5S7 


0.752 


0.722 


0.612 


0.589 


0.536 


0.472 


E th /s [GeV] 


19.25 


10.9 


8.08 


8.21 


7.19 


6.80 



Table 9. Top section: SHM yields of baryon b and at central rapidity db/dy, and strangeness s and ds/dy at RHIC, left 
for the total system, right for the central rapidity region. Next, we give strangeness per entropy s/S (for central rapidity: 
ds/dS), strangeness per baryon s/b, the energy cost to make strangeness pair snn / (2ds / db) , thermal energy per baryon at 
hadronization dE th /db, fraction of initial collision energy in thermal degrees of freedom, (2i? th /6)/ x /SNN- All 62.4 GeV results, 
and the 130 GeV 4n results, are, as discussed in text, result of assumptions, and/or interpolation of yields, and/or parameters, 
and hence are a prediction. 



v^nn [GeV] 


62.4 


130 


200 


62.4 


130 


200 








dN/dy\ y -_ 




b, db/dy 


350 


350 


350 


33.5 


18.5 


14.8 


s, ds/dy 


560 


755 


726 


120.4 


136.7 


123 


WOs/S, ds/dS 


3.17 


2.43 


2.66 


3.30 


3.39 


3.13 


s/b, ds/db 


1.60 


2.16 


2.07 


3.35 


6.87 


8.29 


^s^/(2ds/db) [GeV] 


19.5 


30.1 


48.3 


9.31 


9.46 


12.06 


dE th /db [GeV] 


7.18 


9.52 


11.24 


14.8 


30.4 


38.2 


(2dE th /db)/^/si^ 


0.230 


0.146 


0.112 


0.474 


0.467 


0.382 


dE th /ds [GeV] 


4.49 


4.41 


5.41 


4.42 


4.42 


4.60 



As the bottom right of table 9 indicates, the fraction 
of energy stopped in the central rapidity region at RHIC, 
(2dE t h/db) / ^/snn is rather large, it is estimated to be 58% 
at y / snn = 62.4 GeV decreasing to 36% at top RHIC en- 
ergy. The energy cost to make strangeness extrapolates 
well from the SPS level, connecting smoothly, see the bot- 
tom panel in figure 7, for both total yield and central 
rapidity yield. We note, in passing, that only a small frac- 
tion, 10%, of the total energy is thermalized at the top 
RHIC energy considering the total fireball. 90% is evi- 
dently the energy of the collective flow, predominantly in 
the longitudinal direction. 

The expectation of ever rising strangeness yield with 
y / SNN are not disappointed in Fig. 8, but the rapid smooth 
rise is surprising. One finds such a result in a nearly model 
independent analysis adding up the s carrier particles, 
which are mostly directly measured. A more precise study 
which adds up strangeness in the particles produced ac- 
cording to the SHM as seen in tables 8 and 9 is shown in 
figure 8 — there are non-negligible contributions of unob- 
served hidden strangeness, in particular in the r\ hadron 
(40% ss content). We have scaled the strangeness yield to 
the 7% centrality with iV w = 349 for the total yields. For 



the central rapidity, we present results for the 5% central- 
ity. 



4 Discussion and Interpretation 

4.1 The K+/tt+ horn 

One can wonder how, in qualitative terms, can a param- 
eter 7q, which controls the light quark yield, help explain 
the horn structure seen in figure 2. We observe that this 
horn structure in the K+/-7T+ ratio traces the final state 
valance quark ratio s/d, and in language of quark phase 
space occupancies 7$ and fugacities Aj, we have: 

7T + d \\dj 7d V \) lq 

(11) 

In chemical equilibrium models j s /"fq = 1, and the 
horn effect must arise solely from the variation in the ratio 
\ S /Xq and the change in temperature T. The isospin fac- 
tor A/3 is insignificant in this consideration. For the inter- 
esting range of freeze-out temperature, F(T) is a smooth 
function of T. Normally, one expects that T increases with 
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Fig. 7. The specific strangeness yield as function of reac- 
tion energy -^snn- Top strangeness per baryon s/b, middle 
strangeness per entropy s/S and bottom E t h/s thermal en- 
ergy content per strangeness. Solid squares correspond to 
the triangles on right are for the rapidity density yield dN/dy 
at RHIC. The total yield results are connected by solid line to 
guide the eye, and the central rapidity results (red) are con- 
nected by dashed line. 



collision energy, hence we expect an monotonic increase in 
the K + /7r + ratio, not considering the quark chemistry. 

As collision energy is increased, increased hadron yield 
leads to a decreasing X q = e MB / 3T . We recall the smooth 
decrease of /xb with reaction energy seen in bottom panel 
in figure 5. The two chemical fugacities X s and X q are cou- 
pled by the condition that the strangeness is conserved. 
The chemical potential effect is suggesting a smooth in- 
crease in the K + /ir + ratio. With considerable effort, one 
can arrange the chemical equilibrium fits to bend over 
to a flat behavior at ^/s^ N as the dotted line in figure 
2 shows. It is quasi impossible to generate a horn with 
chemical equilibrium model. 

Consideration of chemical non-equilibrium allows us to 
consider an energy dependent ratio 7 s /7 9 , which as seen 
in Eq. (11) is a multiplicative factor in the horn struc- 
ture. The fit produces a horn like behavior of j s /jq at 
^/s^ N , seen in figure 5. As function of energy, many other 
particle yields must remain relatively smooth, with a few 
exceptions seen in figure 3. We see that the description 
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Fig. 8. The strangeness s (= s) content resulting from the 
SHM fit, as function of reaction energy. The total yield re- 
sults, solid squares (blue), are scaled with number of wounded 
nucleons to a fireball formed in 7% central Pb-Pb reactions 
(TVw = 349). The triangles, on right, are for the rapidity den- 
sity yield ds/dy at RHIC. 

of the horn structure is possible, as there are effectively 
three function of ^/snn which help to create it, T, X q /X s 
and 7s /7g, but it is in no way assured that the right horn 
arises, seen the behavior with energy of the other particle 
yields. 

Indeed, only the full chemical non-equilibrium model 
in which the two phase space occupancies, 7 S and j q , 
vary independently, does a good job as is seen compar- 
ing the solid with dashed and dotted lines in figure 2. 
Seen the horn-like structure of all these lines obtained re- 
laxing strangeness conservation constraint we realize that 
it is not the increased number of parameters, but the fact 
that particle production follows the SHM with chemical 
non-equilibrium which allows the non-equilibrium model 
to succeed. 



4.2 The K+/7T+ horn as function of reaction volume 

The rather sudden changes in freeze-out parameters j q 
and T appears to be a universal behavior. We established 
it here as function of energy, and in earlier work as func- 
tion of the reaction volume (i.e., participant number A), 
see figure 1 in [3]. In both cases, the chemical freeze-out 
temperature is higher below a threshold, as expressed ei- 
ther by low energy or participant number. The most dras- 
tic change is that 7 9 jumps up from a value at, or below 
0.5, to 1.6 as either the energy or volume threshold is 
crossed. The volume threshold is, however, not as sharp 
as the reaction energy threshold. The large system limit is 
achieved for A > 25, with a smooth transition beginning 
at A > 6, as can be seen in figure 4 in [3]. 

Seeing this remark, one immediately wonders if the 
K + /7r + horn is present in the impact parameter kaon and 
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Fig. 9. -\/K+K~/7r+7r _ as function of participant number A 
varying with reaction centrality, PHENIX data [62] . 



pion data and the answer is no. Actually, this is not sur- 
prising: since both ir + and K + originate, in our study, at 
the level of about 50% in directly thermally produced par- 
ticles the ratio K + /ir + is a measure of the horn structure 
is due to a rise in density of strangeness s at hadroniza- 
tion, outpaced by the rise in the d density above ^/ S NN' 
whatever the mechanism in terms of statistical parame- 
ters that implements this. However, when considering the 
impact parameter dependence at v /snn = 200 GeV, the 
rise in strangeness has yet to occur, as in the small volume 
there has insufficient life span to produce strangeness. In 
this situation we do not expect that the horn is present as 
function of A. 

One can see the delayed production of strangeness as 
function of impact parameter directly in the PHENIX 
impact parameter data [62], without need for a detailed 
theoretical analysis. Consider the ratio shown in figure 
9: K/n = *\/K+K~/7r+7r - . This particular product-ratio 
of particles is nearly independent of chemical potentials 
Hb, Ms an( i the volume V since it comprises ratio of prod- 
ucts of particles and antiparticles. The rise seen in fig- 
ure 9 is evidence for an additional strangeness production 
mechanism turning on at about A ~ 20. In figure 9 we 
do not show a common systematic error, thus the normal- 
ization scale of the figure could undergo a revision. This 
cannot change the insight that the additional strangeness 
above and beyond the first collision content is produced for 
A > 20, enhancing the global yield by 50% or more. More- 
over, we see that the rise is gradual as can be expected in 
kinetic theory models of strangeness production [82,84], 
and there is at the maximum centrality no evidence as yet 
of strangeness yield saturation. 

The entropy content of the small system A < 20 is such 
that strangeness per entropy is at the level of s/S ~ 0.02, 
and both entropy and strangeness rise with centrality of 
the reaction at -^/snn = 200 GeV. However, unlike the 
energy dependence, the ratio s/S rises modestly, strang- 
eness does not outpace entropy rise by more than 20%. 
This is in agreement with expectation, since the threshold 
of strangeness mass is not fully relevant at the top RHIC 



reaction energy, and thus we are seeing the properties of 
a deconfined initial state in which strange quark is effec- 
tively massless. Instead, it is the lifespan of the system 
that matters, as noted above. 

There is very little observed dependence of ratios of 
hadron resonances with the ground state yields, such as 
K* / K. This implies and agrees in quantitative way with 
the tacit assumption inherent in the above discussion, and 
the result of a more detailed analysis [3] , that there is no 
T dependence of the freeze-out conditions for A > 20. For 
this reason for A > 20 ratios of all hadrons which do not 
involve a difference in strangeness content, do not vary 
with centrality. 

We further note that there is little change in chemical 
potentials with centrality for A > 20, indicating that the 
stopping of baryons is not a result of multiple scattering, 
but is due to phase conditions of matter. Comparing other 
properties of matter, we see very much the same behav- 
ior as function of impact parameter and reaction energy: 
in particular, we note the step up in pressure, in energy 
density, and in entropy density at the impact parameter 
threshold [3]. 



4.3 Chemical equilibrium or non-equilibrium? 

An important questions discussed in the study of hadron 
yields interpretations is if chemical equilibrium or non- 
equilibrium prevails in the hadronization process. There 
are workers who strongly defend the chemical equilibrium 
SHM [85]. Let us look again at the survey of fit quality 
results seen in figure 1. We note that for j q — 1 (but 
7 S =/= 1) at each energy there seems to exist a reasonable 
fit with 0.5 < % 2 /d.o.f. < 1.5 for the data sample con- 
sidered, which suggests that at each reaction energy with 
7 9 = 1 a reasonable and widely accepted physical descrip- 
tion of the experimental data emerges. This result is there- 
fore claimed in studies that focus on the hadron yields at 
each energy separately. What works poorly in SHM used 
with 7 9 = 1 and even worse with 7 g = 1, 7 S = 1 is the 
energy dependence of particle ratios, with the most promi- 
nent present day example being the horn structure in the 
K + /7r + yield ratio. Seen from this perspective, it is the en- 
ergy dependent particle yield that requires the inclusion 
in the necessary set of parameters a varying value "f q ^ 1. 

Another important question directly related to the is- 
sue of chemical equilibrium and non-equilibrium is how 
the fitted results for T(/xb), the 'hadronization curve' re- 
late to the phase boundary between deconfined primary 
phase and the hadron phase. Clearly, the result of the fit 
are greatly dependent on the assumption about chemical 
condition with the equilibrium fit claiming a hadroniza- 
tion at RHIC at T = 175 MeV. 

The rapidly decreasing freeze-out temperature T as 
•y/sNN decreases, and which is certainly inconsistent with 
the rather flat phase transition boundary at moderate 
chemical potentials is explained by suggesting that the 
hadronization may be related to a particular values of en- 
ergy per particle content, of the magnitude 1 GeV [85]. 
However, this condition, though not rooted in any known 
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basic physical principle, is also obtained in some dynami- 
cal studies, see, e.g., Ref. [86,87]. We note that the chem- 
ical equilibrium hypothesis fails to explain the hadroniza- 
tion conditions expected as function of T and /is , or equiv- 
alently, as function of ^/snn- 

In summary, the interpretation of hadron production 
in terms of chemical equilibrium SHM disagrees, in quan- 
titative manner, both with the reaction energy dependent 
particle yields (such as the K + /ir + horn) and the reaction 
energy dependent shape of hadronization boundary. 



4.4 Hadronization boundary in heavy ion collisions 

We believe that the hadronization boundary, in the T- 
fiB plane, is the result of a complex interplay between 
the dynamics of heavy ion reaction and the properties 
of both phases of matter, the inside of the fireball, and 
the hadron phase we observe. Even disregarding compli- 
cations related to the rapid expansion of the dense matter 
fireball, the presence of chemical non-equilibrium particle 
distributions introduces significant freedom into the shape 
and location of the T(^b) transition region. 

Recall, first, that available lattice results apply to a 
system in the thermodynamic limit with j q = j s = 1, 
for both quark and confined hadron phases. The typical 
boundary between the QGP and hadron phases is dis- 
cussed in Ref. [77], and is dependent on chemical proper- 
ties of QGP. Typically, one considers the dependence on 
chemical potentials, and in particular on hb, however, a 
significant change in the phase boundary location is to 
be expected when j q and j s + 1 1. To understand this 
important remark, consider the two other known cases 
7q = 1, 7 S = corresponding to 2 flavors, and j q = j s = 
corresponding to pure gauge. There is a significant change 
in T(/ib = 0), which increases with decreasing 7$. 

Moreover, not only the location but also the nature of 
the phase boundary can be modified by variation of 7$. We 
recall that for the 2+1 flavor case, there is possibly a crit- 
ical point at finite baryochemical potential with — 350 
McV [77,88]. However, for the case of 3 massless flavors 
there can be a 1st order transition at all /ib [89,90]. Con- 
sidering a classical particle system, one easily sees that an 
over-saturated phase space, e.g., with -f q — 1.6, 7 S > "f q 
for the purpose of the study of the phase transition acts as 
being equivalent to a system with 3.2 light quarks and 1.6 
massive (strange) quarks present in the confined hadron 
phase. 

Even though one should be keenly aware that over- 
saturation of the phase space is not the same as additional 
degeneracy due to true degrees of freedom, the similarity 
of resulting effect must be considerable. We know that 
with increasing /j,b, the increased quark density creates 
the environment in which the phase cross-over becomes 
a phase transition. The influence of 7^ p cannot be dif- 
ferent. Considering that 7^ p enhances both quark and 
antiquark number, it should be more effective compared to 
/is in its facilitation of a phase transition, and reduction 
of the temperature of the phase boundary for 7^ p > 1 • 



We therefore can expect that, for a chemically over- 
saturated system, there is also an effective increase in the 
number of degrees of freedom. Looking at the structure 
of the quark-hadron transformation this increase in the 
number of available effective degrees of freedom occurs in 
a physical system which is almost, but not quite, able to 
undergo a 1st order phase transition. Considering here also 
the sudden nature of the fireball breakup seen in several 
observables [1], we conjecture that the hadronizing fireball 
leading to 7 S > "f q — 1.6 passes a true phase boundary 
corresponding to a 1st order phase transition condition at 
small /ib- Because of the changed count in the degrees of 
freedom, we expect that the phase transition temperature 
is at the same time decreased to below the cross-over value 
for chemical equilibrium case of 2+1 flavors near T = 
162 ± 3 MeV. 

It seems to us that it would be very interesting to 
determine, in as more rigorous way for the case of the 
2+1 flavor lattice QCD at ^b = for which values, if 
any, of 7$ the system undergoes a phase transition of 1st 
order. Lattice QCD methods employed to obtain results 
at finite /ib, e.g., the power expansion [88,91], should also 
allow to study the case of y(i 7 = Tlnji > 0, and near to 
/Lt 7 =, i.e. 7^ = 1. We see the actual difficulty in the need 
to simulate different values of /i 7 in the two phases, such 
that the quark pair content is preserved across the phase 
boundary. 

The dynamical, and theoretically less spectacular, ef- 
fect capable to shift the location in temperature of the ex- 
pected phase boundary, is due to the expansion dynamics 
of the fireball. The analysis of the RHIC results suggests 
that the collective flow occurs at parton level [92] . Collec- 
tive flow of color (partons) is like a wind capable to push 
out the color non transparent 'true' vacuum [93], adding 
to thermal pressure the dynamical component, for a finite 
expanding system this would lead to supercooling [94]. 
This dynamical effect will push the hadronization condi- 
tion to lower local freeze-out T at high ^/snn, thus flatten- 
ing the boundary between the phases as function of /ib ■ In 
the context of results we have obtained, it is the smooth- 
ness of the ratio P/e obtained at hadronization which sup- 
ports the possible relevance of dynamic phase boundary 
displacement. This behavior suggests a smoothly chang- 
ing dynamical break up condition, potentially related to 
(hydrodynamic) flow. 

4.5 Our hadronization boundary and its interpretation 

The above two effects, the change in the location of the 
static phase boundary in presence of chemical non-equili- 
brium and the dynamics of collective matter flow toward 
breakup condition, are both non-negligible but hard to 
evaluate quantitatively. We believe that they can explain 
why the chemical freeze-out conditions T and hb are as 
presented in figure 5. Of particular relevance is the low 
value of T at high reaction energy, and relatively high 
value of T at low reaction energy, just opposite what one 
finds to be result of SHM analysis when chemical equilib- 
rium is assumed. 
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Fig. 10. T—fiB plane with points obtained in the SHM fit. See 
text for discussion. 



Our chemical freeze-out conditions are better shown in 
the T-/xb plane, see figure 10. Considering results shown 
in figure 5, we are able to assign to each point in the T-/iB 
plane the associate value of v /snn. The RHIC dN/dy re- 
sults are to outer left. They are followed by RHIC and SPS 
iV 4w results. The dip corresponds to the 30 and 40 AGeV 
SPS results. The top right is the lowest 20 ,4GeV SPS 
and top 11.6 ^4GeV AGS energy range. We see that the 
chemical freeze-out temperature T rises for the two lowest 
reaction energies 11.6 and 20 A GeV to near the Hage- 
dorn temperature, T = 160 MeV. Such phase structure is 
discussed e.g. in the context of chiral quark pairing [12]. 

The size of error bars, in figure 10, is output of the fit 
process, and when it is rather large, it implies that the 
resulting x 2 profile was relatively flat, or/and that there 
were two neighboring god fit minima. To guide the eye, we 
have added two lines connecting the fit results. As seen in 
figure 10, at //b = 0, we find that hadronization occurs 
at T = 140, decreasing to T = 120 MeV at fi B = 400 
MeV. Along this line 7 9 > 1.6. As argued in the previous 
subsection 4.4, this line could be a true 1st order phase 
boundary between quark matter and an over-saturated 
hadron phase. 

Two different interpretations come to mind when we 
attempt to understand the other branch in figure 10, the 
rise from T ~ 120 to 160 for /is > 400, which is ac- 
companied by a rather low hadron side phase space occu- 
pancy. Most 'natural' is to presume that the dissolution 
of color bonds did not occur in heavy ion collisions be- 
low 30 A GeV, we are dealing with 'conventional' hadron 
matter. The under-saturation occurs since there was no 
time to make hadrons, i.e., chemical equilibration was not 
achieved in the colliding hadron system before it breaks 
apart. 

The fireball break-up at a higher temperature is a 
consistent freeze-out scenario for under-saturated hadron 
phase space considering the kinetic scattering freeze-out 
condition. Given the greatly reduced particle density (oc 
7™, n = 3, 2, 1) a high T freeze-out for 7,; < 0.5 is consis- 



tent. The nucleon density scales with 7^ and pion density 
with 7^. Therefore the meson-baryon and meson-meson 
scattering length scales as L oc I/7 5 , and L oc I/7 4 , re- 
spectively. This implies that, as the system expands, it is 
able to freeze out early at a higher T. 

On the other hand, the volume size we found, see ta- 
ble 1, is significantly larger at low reaction energies. This 
implies that a scenario with pure hadron matter present is 
subject to a quite different expansion history. This signals 
that a standard picture of a conventional hadron matter 
formation at reaction energies below the transition point 
at 6.26 GeV < y/s^ < 7.61 GeV may not be the valid 
explanation of results of our analysis. Namely, if the re- 
action history since first contact is different for the two 
reaction energy ranges, one would expect that the sys- 
tematics of the final state entropy production, strangen- 
ess production, and strange antibaryon production has a 
visible break at the critical point. What we have found 
is, instead, that these quantities show a rather smooth 
uninterrupted rise with reaction energy. 

This means that the the initial conditions reached in 
the reaction where, e.g., entropy and strangeness are pro- 
duced, is not undergoing a sudden change. The change oc- 
curs at the end near to the hadronization of matter. For 
this reason, we see a change in particle yields (the horn), 
statistical parameters jump, and the physical conditions 
at hadronization jump even more. The yields of quantities 
which are driven by physics of the initial dense matter for- 
mation, e.g., entropy content, strangeness content, change 
smoothly with heavy ion reaction energy in the domain 
we explored. We are furthermore swayed away from the 
picture of the hadronic gas being the form of matter at 
breakup below ^/^nn ty the strange antibaryon produc- 
tion systematics we discussed in figure 4. 

We are searching thus for an explanation in terms of a 
new phase of matter being involved in the hadronization 
process, but clearly this cannot be the semi-perturbative 
quark-gluon plasma state. The conceivable explanation 
of the fit result below 30 A GeV is presence, at the high 
baryon density arising at large [1b , of a constituent quark 
plasma [95]. Even if the perturbative QCD quark phase 
is reached at high temperature, in expansion-cooling the 
system encounters the valon (word derived from 'valance' 
quark) phase in which the color quark bonds are broken, 
but chiral symmetry restoration is not completed, with 
quarks of mass m U: d — 340 MeV and m s ~ 500 MeV 
being the only active degrees of freedom. This scenario is 
not inconsistent with the finding on the lattice, that, for 
/ie and in chemical equilibrium, the chiral symmetry 
restoration coincides with deconfinement transition. 

In a valon matter, even assuming chemical equilibrium, 
the number of quark pairs at temperature near to T = 160 
MeV would be rather small, given the high constituent 
quark mass. In breakup of this system, a relatively small 
7^ G is achieved. Furthermore, since the mass of these con- 
stituent quarks is greater than that of the pion, the phase 
transformation between hadron and valon matter occurs 
at relatively large T. To see this, recall that the pion 
with its low mass produces greater pressure than valons 
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and thus is pushing the transition boundary to higher T. 
Strangeness, and importantly the entropy content in this 
phase arise due to prior initial state perturbative QGP 
phase and hence such a valon system must be larger in 
volume at the point of hadronization. 

It is also conceivable that a hadron fireball evolving 
from the beginning and fully in the valon phase would 
maintain much of the continuity we saw in hadronic ob- 
servables. For example, u and c?-valon-quark scattering can 
produce strange valon-quark pairs, and these give rise in 
hadronization to the abundances of strange antibaryons, 
as expected in the deconfined phase. What speaks for this 
option is the rather sudden change in the thermal energy 
content per strange quark pair produced, which is seen 
in bottom of figure 7, indicating appearance of a new en- 
ergy efficient mechanism of strangeness production above 

/ c cr 
V^NN- 

One may wonder how our findings compare to earlier 
studies of the phase boundary, both in statistical mod- 
els [96], and microscopic models (see Ref. [87] and ref- 
erences therein). In the microscopic models one accom- 
plishes a better understanding of the approach to ther- 
mal and chemical equilibrium of the degrees of freedom 
employed. A continuous phase boundary is, here, a di- 
rect outcome of the assumption made about the degrees 
of freedom present. Our analysis, which does not rely on 
such assumptions, is thus less model dependent and al- 
lows for presence of degrees of freedom with unexpected 
properties. On the other hand, we also firmly believe, as 
is shown in figure 10 that there is a smooth phase bound- 
ary, with T dropping with increasing /xb ■ What our study 
has uncovered is the possible presence of another phase 
boundary for /ib > 350 MeV at higher T. It is impor- 
tant for the reader to keep in mind that this finding is not 
in conflict with theoretical chemical equilibrium results 
which focus on the other, conventional, phase branch and 
address physics of phase transformation occurring in the 
early Universe. 

Moreover, a recent study of the low energy AGS pion 
production data [97] found that the thermal freeze-out 
temperature at reaction energy of 8 GeV is at T t h = 140 
MeV [98]. Thus there is also consistency of our present 
analysis with the shape of pion transverse mass spectra 
and the high chemical freeze-out temperature we find at 
AGS. Another interesting finding was the medium mass 
modification which allowed to describe pion decay spectra. 

4.6 Final remarks 

In summary, we have performed a complete analysis of 
the energy dependence of hadron production in heavy ion 
collisions, spanning the range beginning at the top AGS 
energy, to the top RHIC energy. We have made extensive 
predictions about particle production in the entire energy 
range. These results are useful in several respects. For ex- 
ample, we have shown that the best energy to search for 
the elusive pentaquarks would be at SPS at 30-40 A GeV, 
where we find that the total yield of 6> + (1540) is already 
fully developed. Thus there is a maximum in the ratio 



+ /K s ~ 0.2 at 30 A GeV. Of course, this finding pre- 
supposes the existence of the exotic state. 

We have furthermore presented hadron yields impor- 
tant in the understanding of dilepton spectra, such as p, rj 
and u>. The relative meson resonance yields we find do not 
follow the pp systematics and vary as function of energy. 
One thus can test the hadronization picture here presented 
in study of resonance production. This observation was re- 
cently exploited in a systematic fashion [99] . 

We have shown that the threshold in energy which 
generates a horn in the K + /tt + yield ratio can be associ- 
ated with the chemical freeze-out shifting rather rapidly 
toward condition of greatly increased hadronization den- 
sities. This transition separates the high entropy density 
phase at high heavy ion reaction energy from a low en- 
tropy density phase. This behavior parallels the findings 
for impact parameter dependence of RHIC results, where 
the low entropy density phase is seen for small reaction 
volumes present at large impact parameters [3]. 

Several observables, including strangeness production, 
show continuity across the energy threshold at 6.26 GeV < 
^/s£j r N < 7.61 GeV, thus, it seems that the critical con- 
ditions expresses a change in the nature of the fireball 
breakup, and to a lesser degree a reaction energy depen- 
dent change in the nature of initial conditions reached in 
the reaction. 

We have discussed, in depth, our findings about the 
hadronization condition T([ib) and have argued that at 
high reaction energies a lst-order phase transition is aris- 
ing in the chemically non-equilibrated hot hadronic matter 
system. Detailed discussion was presented about possible 
changes in phases of hadronic matter as function of reac- 
tion energy and reaction volume. 

This manuscript nucl-th/0504028, on which this work is based, 
has been first web-published in April, 2005. We undertook the 
current revision to correct an entropy yield error which the 
early SHARE release contained. This changes entropy S related 
table entries and figures, the data fits are unaffected. When re- 
doing the fits to obtain entropy, we incorporated the latest 
strange hadron yields of NA49 as available. We thank NA49 
and PHENIX collaborations for valuable comments regard- 
ing the acceptances of weak decays. We thank M. Gazdzicki 
and G. Torrieri for valuable comments. Work supported in 
part by a grant from: the U.S. Department of Energy DE- 
FG02-04ER4131. LPTHE, Univ. Paris 6 et 7 is: Unite mixte 
de Recherche du CNRS, UMR7589. 



References 

1. For theoretical evaluation see: 

Nucl. Phys. A 750, Issue 1, pp 1-171 (March 2005); 
for experimental evaluation see: 

Nucl. Phys. A 757, Issue 1, pp 1-283 (August 2005). 

2. M. Gazdzicki et al. [NA49 Collaboration], J. Phys. G 30, 
S701 (2004). 

3. J. Rafelski, J. Letessier and G. Torrieri, "Centrality de- 
pendence of bulk fireball properties at RHIC," [arXiv:nucl- 
th/0412072], Phys. Rev. C (in press) (2005). 



22 



Jean Letessier, Johann Rafelski: Hadron Production and Phase Changes in Relativistic Heavy Ion Collisions 



4. J. Letessier and J. Rafelski, Phys. Rev. C 59, 947 (1999). 

5. J. Kapusta, B. Muller, J. Rafelski, "Quark-gluon plasma : 
Theoretical foundations, an annotated reprint collection", 
Elsevier (Amsterdam 2003), and references therein. 

6. J. Rafelski and B. Muller, Phys. Rev. Lett. 48, 1066 (1982) 
[Erratum-ibid. 56, 2334 (1986)]. 

7. N. K. Glendenning and J. Rafelski, Phys. Rev. C 31, 823 

(1985) ; note that in Fig. 2 the pion yield from resonance 
decays is not included, this dilutes the expected yields by 
about factor two, flattening the curve. 

8. J. Letessier, A. Tounsi, U. Heinz, J. Sollfrank and J. Rafel- 
ski, Phys. Rev. Lett. 70, 3530 (1993), and 

Phys. Rev. D 51, 3408 (1995). 

9. P. Koch, B. Muller and J. Rafelski, Phys. Rept. 142, 167 

(1986) . 

10. P. Koch and J. Rafelski, Nucl. Phys. A 444, 678 (1985). 

11. J. Rafelski and M. Danos, Perspectives In High-Energy Nu- 
clear Collisions, NBSIR-83-2725 and GSI-83-6, see section 
5, KEK scan available at 

http://ccdb3fs.kek.jp/cgi-bin/img_index7200031578; 
update: Lecture Notes in Physics, 231, 361 (1985) 
( Springer- Verlag, Berlin); 

M. Danos and J. Rafelski, Heavy Ion Phys. 14, 97 (2001) 
(M. Danos memorial volume). 

12. M. Kitazawa, T. Koide, T. Kunihiro and Y. Nemoto, Phys. 
Rev. D 65, 091504(R) (2002). 

13. M. Gazdzicki, Acta Phys. Polon. B 34 (2003) 5771; 

M. Gazdzicki and M. I. Gorenstein, Acta Phys. Polon. B 
30, 2705 (1999). 

14. G. Torrieri, W. Broniowski, W. Florkowski, J. Letessier 
and J. Rafelski, Comp. Phys. Com. 167, 229 (2005), and 
references therein by the authors; and 

G. Torrieri, S. Jeon, J. Letessier and J. Rafelski, Comput. 
Phys. Commun. 175 (2006) 635, and see further updates 
at: www.physics.arizona.edu/~torrieri/SHARE/share.html 

15. G. Torrieri and J. Rafelski, J. Phys. G 30, S557 (2004). 

16. J. Letessier and J. Rafelski, Int. J. Mod. Phys. E 9, 107, 
(2000). 

17. P. Koch, J. Rafelski and W. Greiner, Phys. Lett. B 123, 
151 (1983). 

18. I. G. Bearden et al. [BRAHMS Collaboration], Phys. Rev. 
Lett. 90, 102301 (2003). 

19. F. Becattini, J. Cleymans, A. Keranen, E. Suhonen and 
K. Redlich, Phys. Rev. C 64, 024901 (2001). 

20. V. Friese [NA49 Collaboration], Nucl. Phys. A 698, 487 
(2002). 

21. J. Rafelski, J. Letessier and G. Torrieri, Phys. Rev. C 64, 
054907 (2001) [Erratum-ibid. C 65, 069902 (2002)]. 

22. C. Blume [NA49 Collaboration], J. Phys. G 31, S685 
(2005). 

23. J. I. Kapusta and S. M. H. Wong, Phys. Rev. Lett. 86, 
4251 (2001). 

24. J. Letessier, J. Rafelski and G. Torrieri, "Deconfinement 
energy threshold: Analysis of hadron yields at 11.6-A- 
GeV," [arXiv:nucl-th/0411047]. 

25. The results of NA49 experiment used here have been 
provded upon our request by M. Gazdzicki and B. Lung- 
witz, (private communication, compilation of September 
2006), which documentation includs explanation of the 
treatment of weak decays. Among sources of these results 
we note in particular the follwoing works [2,20,22,26-51]. 
The abstraction of these 28 works with in part not entirely 
consistent result tables (due to varying WD corrections, 



and some time dependence of preliminary results) into a 
11 page table set appears to us to be solely possible by 
these key members of the NA49 collaboration. We thank 
Drs. M. Gazdzicki and B. Lungwitz, for their kind help. 

26. V. Friese et al. [The NA49 Collaboration], J. Phys. G 30, 
S119 (2004). 

27. C. Alt et al., Phys. Rev. C 73 (2006) 044910. 

28. A. Richard [NA49 Collaboration], J. Phys. G 31 (2005) 
S155. 

29. C. Blume [NA49 Collaboration], J. Phys. G 31 (2005) S685. 

30. C. Alt et al. [The NA49 Collaboration], J. Phys. G 30, 
S119 (2004). 

31. P. Dinkelaker [NA49 Collaboration], J. Phys. G 31 (2005) 
S1131. 

32. S. V. Afanasiev et al. [The NA49 Collaboration], Phys. 
Rev. C 66 (2002) 054902. 

33. T. Anticic et al., Phys. Rev. C 69 (2004) 024902. 

34. S. V. Afanasiev et al., Nucl. Phys. A 715 (2003) 161. 

35. T. Anticic et al. [NA49 Collaboration], Phys. Rev. Lett. 
93, 022302 (2004). 

36. C. Meurer [NA49 Collaboration], J. Phys. G 30, S1325 

(2004) . 

37. M. Mitrovski [NA49 Collaboration], Acta Phys.Hung. A24 

(2005) 157.. 

38. C. Alt et al. [NA49 Collaboration], Phys. Rev. Lett. 94 

(2005) 192301. 

39. J. Bachler et al. [NA49 Collaboration], Nucl. Phys. A 661 
(1999) 45. 

40. C. Alt et al. [NA49 Collaboration], Eur. Phys. J. C 45 

(2006) 343. 

41. S. V. Afanasev et al. [NA49 Collaboration], J. Phys. G 27 
(2001) 367. 

42. S. V. Afanasev et al. [NA49 Collaboration], Phys. Lett. B 
491 (2000) 59. 

43. C. Alt et al. [NA49 Collaboration], "Pion and kaon pro- 
duction in central Pb+Pb collisions at 20A and 30A GeV: 
Evidence for the onset of deconfinement," arXiv:0710.0118 
[nucl-ex] . 

44. T. Susa [NA49 Collaboration], Nucl. Phys. A 698 (2002) 
491. 

45. C. Alt et al. [NA49 Collaboration], Phys. Rev. Lett. 94 
(2005) 052301. 

46. I. Kraus [NA49 Collaboration], J. Phys. G 31 (2005) S147. 

47. A. Mischke et al., Nucl. Phys. A 715, 453 (2003). 

48. G. I. Veres [NA49 Collaboration], Nucl. Phys. A 661 (1999) 
383. 

49. V. Friese [NA49 Collaboration], Nucl. Phys. A 698 (2002) 
487. 

50. S. V. Afanasiev et al. [NA49 Collaboration], Phys. Lett. B 
538 (2002) 275. 

51. M. Mitrovski, J.Phys. G32 (2006) S43. 

52. J. Rafelski and J. Letessier, Acta Phys. Polon. B 34, 5791 
(2003); and J. Phys. G 30 (2004) SI. 

53. B. B. Back et al. [E917 Collaboration], Phys. Rev. C 69, 
054901 (2004). 

54. J. Cleymans, H. Oeschler, K. Redlich and S. Wheaton 
Phys. Lett. B 615, 50 (2005). 

55. J. Cleymans and K. Redlich, Phys. Rev. C 60, 054908 
(1999). 

56. K. Adcox et al. [PHENIX Collaboration], Phys. Rev. C 
69, 024904 (2004). 

57. C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 89, 
092301 (2002). 



Jean Letessier, Johann Rafelski: Hadron Production and Phase Changes in Relativistic Heavy Ion Collisions 



23 



58. K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 
89, 092302 (2002). 

59. J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 
92, 182301 (2004). 

60. C. Adler et al. [STAR Collaboration], Phys. Rev. C 66, 
061901 (2002). 

61. C. Adler et al, Phys. Rev. C 65, 041901 (2002). 

62. S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. C 
69, 034909 (2004). 

63. H. B. Zhang [STAR Collaboration], "Delta, K* and rho 
resonance production and their probing of freeze-out dy- 
namics at RHIC," [arXiv:nucl-ex/0403010]; 

J. Adams [STAR Collaboration], Phys. Rev. C 71 064902 
(2005). 

64. C. Markert [STAR Collaboration], J. Phys. G 30, S1313 
(2004). 

65. J. Adams et al. [STAR Collaboration], Phys. Lett. B 612, 
181 (2005). 

66. D. Elia [NA57 Collaboration], J. Phys. G 31, S135 (2005). 

67. I. G. Bearden et al. [BRAHMS Collaboration], Phys. Rev. 
Lett. 94, 162301 (2005). 

68. B. B. Back et al, Phys. Rev. Lett. 91, 052303 (2003). 

69. I. G. Bearden et al. [BRAHMS Collaboration], Phys. Rev. 
Lett. 93, 102301 (2004). 

70. J. Rafelski, Phys. Rept. 88, 331 (1982). 

71. G. E. Bruno [NA57 Collaboration], J. Phys. G 30, S717 
(2004). 

72. C. Alt et al. [NA49 Collaboration], Phys. Rev. Lett. 92, 
042003 (2004). 

73. A. R. Dzierba, D. Krop, M. Swat, S. Teige and 
A. P. Szczepaniak, Phys. Rev. D 69 (2004) 051901 
[arXiv:hep-ph/0311125]. 

74. J. Letessier, G. Torrieri, S. Steinke and J. Rafelski, Phys. 
Rev. C 68, 061901 (R) (2003). 

75. Y. Aoki, G. Endrodi, Z. Fodor, S. D. Katz and K. K. Szabo, 
Nature 443 (2006) 675 [arXiv:hep-lat/0611014]. 

76. Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. Lett. 
B 643 (2006) 46 [arXiv:hep-lat/0609068]. 

77. Z. Fodor and S. D. Katz, JHEP 0404, 050 (2004), and 
references therein. 

78. C. R. Allton et al., Phys. Rev. D 66, 074507 (2002). 

79. I. Kuznetsova and J. Rafelski, Eur. Phys. J. C 51 (2007) 
113 [arXiv:hep-ph/0607203]. 

80. J. Rafelski and M. Danos, Phys. Lett. B 97, 279 (1980). 

81. J. Rafelski and J. Letessier, J. Phys. G 28, 1819 (2002). 

82. J. Letessier, A. Tounsi and J. Rafelski, Phys. Lett. B 389 
(1996) 586. 

83. J. Rafelski and J. Letessier, Phys. Lett. B 469, 12 (1999). 

84. J. Letessier and J. Rafelski, Phys. Rev. C 75 (2007) 014905 
[arXiv:nucl-th/0602047]. 

85. P. Braun-Munzinger, K. Redlich and J. Stachel, in Quark 
Gluon Plasma 3, eds. R.C. Hwa and Xin-Nian Wang, 
World Scientific Publishing, and references therein. 

86. M. Bleicher and J. Aichelin, Phys. Lett. B 530, 81 (2002). 

87. E. L. Bratkovskaya et al, Phys. Rev. C 69, 054907 (2004). 

88. C. R. Allton, M. Doring, S. Ejiri, S.J. Hands, O. Kacz- 
marek, F. Karsch, E. Laermann, and K. Redlich, Phys. 
Rev. D 71, 054508 (2005). 

89. A. Peikert, F. Karsch, E. Laermann and B. Sturm, Nucl. 
Phys. Proc. Suppl. 73, 468 (1999). 

90. C. Bernard et al. [MILC Collaboration], Phys. Rev. D 71, 
034504 (2005). 



91. C. R. Allton, S. Ejiri, S. J. Hands, O. Kaczmarek, 
F. Karsch, E. Laermann and C. Schmidt, Phys. Rev. D 
68, 014507 (2003). 

92. H. Z. Huang and J. Rafelski, in AIP Conf. Proc. 756, 210 
(2005). 

93. T. Csorgo and J. Zimanyi, Heavy Ion Phys. 17, 281 (2003). 

94. J. Rafelski and J. Letessier, Phys. Rev. Lett. 85, 4695 
(2000). 

95. I. I. Roizen, E. L. Feinberg and O. D. Chernavskaya, Phys. 
Usp. 47, 427 (2004) [Usp. Fiz. Nauk 47, 427 (2004)]. 

96. F. Becattini, M. Gazdzicki, A. Keranen, J. Manninen and 
R. Stock, Phys. Rev. C 69, 024905 (2004). 

97. J. L. Klay et al. [E-0895 Collaboration], Phys. Rev. C 68, 
054905 (2003). 

98. B. Pin-zhen and J. Rafelski, Eur. Phys. J. A 32, 267 
(2007). 

99. G. Torrieri and J. Rafelski, Phys. Rev. C 75, 024902 (2007) 
[arXiv:nucl-th/0608061]. 



